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On Two-Dimensional Variational Problems 


in Parametric Form 


H. B. JENKINS 


Communicated by R. FINN 


I. Introduction. We shall consider two-dimensional variational problems of 
the type 


(1) 6f £(p)dsdt=0. 


The integral applies to surfaces given in parametric form by functions x(s, #), 
WiSit)pee(St),, and p=(o,7,m) where. o=0(y, 2)/O(s, t), t= Cl(z,x)/0(s, 0), 
w =0(x, y)/0(s,¢). We shall assume that Y£(p) belongs to the Hélder class C?** for 
allp=-0, and that (i) 4 (p) is positive homogeneous of degree 1 ;7.¢., J(Ap) =A Y(p) 
forallA>0;and (ii) Y£(p) is positive definite ; 7. e., there exist positive constants m,, 
my such that m,|p|< 4(p)<m,|p|. For applications we shall assume in addition 
the condition (111) Y£(p) is regular. We define Y(p) to be regular if the surface Y 
defined by Y (p) =1 has everywhere positive Gauss curvature. This definition of reg- 
ularity is equivalent to the condition that the Euler equations of the corresponding 
local non-parametric variational problems are elliptic. Note that conditions (i) and 
(ii) imply that 2'is a closed surface, star-shaped about the origin, since if #(p)=1, 
then p=WN/ Y(N), where N is a unit vector in the direction of p. Our results 
apply to solutions of the Euler equations of (1), hence we shall always use the 
term “‘solution’’, or “‘solution surface’’, in that sense. If N(q) is the unit normal 
at the point q of a solution surface S, p(q)=N(q)/4(N(q)) defines a mapping 
of S into & which is a generalization of the classical spherical image mapping. 
We shall show that the behavior of S in a neighborhood of a point qg is governed 
by the second fundamental form of 2’ at the image point p(q). 

It is our purpose here to study the geometrical relationships between the 
surface 2 and a solution surface. We shall show that the mapping p(q) determines 
a non-trivial linear mapping of the differential forms on 2 into those on a solution 
surface S; and that in the case of regular variational problems this mapping in 
turn gives rise to a canonical mapping of S onto the solution of a non-parametric, 
uniformly elliptic variational problem. 

The principal application is to show that solutions of variational problems 
which satisfy conditions (i), (ii), (iii) behave qualitatively like minimal surfaces, 
in the same sense that solutions of uniformly elliptic equations behave like 
solutions of the Laplace equation. We shall show that the mapping p(q) of a 
solution surface S into XY possesses a property which is a natural generalization 

Arch, Rational Mech. Anal., Vol. 8 13 


182 H. B. JENKINS: 


of the conformality of the spherical mapping of a minimal surface into the unit 
sphere (Theorem 1). It is an immediate consequence of our Besleoss that the 
quasi-linear equations satisfied by solutions of the form z = (x, y) are equations 
of minimal surface type” in the sense of Frnn [/] and that they can be written . 
in the form (A (wu, v)),+(B(u, v)), =0, where u=9,, V=Qy, and A and B are 
bounded functions*. Theorems of Finn [7], [2], [3] may therefore be applied 
to show that such solutions share with non-parametric minimal surfaces the 
properties (a) the gradient of p at a point (x, vg) can be estimated in terms of 
the distance 7 to the boundary anda bound M on the quantity | p(x, vy) —p(%o, Yo) |**; 
(b) there are no solutions with isolated singularities; and (c) the boundary value 
problem can be solved for continuous data on convex curves. Again in the case 
of solutions of the form z=@(x, y), a theorem of the author [4], can be applied 
to show that the gradient of a solution defined in the exterior of a circle goes 
Hélder continuously to a finite limit at infinity. In the minimal surface case, 
the latter property and property (b) were established by Bers [4], and (a) and (c) 
were given by FINN. 

Let S be a simply connected minimal surface whose normal direction at each | 
point makes an angle =f with a fixed direction. Let g be a point of S such that 
the distance along S to the boundary is =d. Let « be the angle between the normal 
at g and the fixed direction. OssERMAN [7] has recently shown that the Gauss 
curvature K(q) of S at g then satisfies an inequality 


1 [4(v?—a?) h(a,7) |? 
(2) IK@| s+ 


7(1—a?)2 : 
where 7 =cot f/2, a=cot «/2, and /(a,7) is a known function. We shall show — 
that the Gauss curvature of a solution of a regular parametric variational problem 
of the type (1) satisfies an inequality 


(3) piece ise 


where «, 8, g, d are as defined above, and C is a constant which depends only on 
B and Y(p). This inequality is qualitatively similar to OSSERMAN’s inequality 
and may be interpreted as implying that the curvature at a point g of a solution 
surface S is small not only when the distance to the boundary is large, but also 
when the spherical image of the point q is close to the boundary of the spherical 
image of S. In the case of solutions of the form z=q(x, y), the inequality (3) 
may be put in the form 


(4) 1K, ¥0)| < ~& 


et ee 

Wo? y2? 

where ¢ is the radius of the largest circle about the point (x, Vo) In which @ is 
defined, and Wy =(1+ q5(%, Yo) +95 (%o. %o))*: A characteristic property of 
solutions of equations of minimal surface type is that the spherical mapping is 
quasi-conformal, and this in turn implies that the ratio of the magnitudes of the 


* If an equation can be written in the form (A (u, v)),+(B(u, v)),=0, then the | 
pair (A, B) is said to be a conservation law of the equation (see LOEWNER [6]). 

** It will be shown in another paper that an estimate of the gradient in terms of 7 
and M can be obtained from the inequality (4), below, by elementary means. 
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principal curvatures is bounded above and below by positive constants. Combining 
the latter property with (4) results in 


(5) +25, +1 < 


where % =,4(%o, Yo)» So =Pxy(Xo» Vo)» 9 =Pyy(%o, Yo), and W, and 7 are as de- 
fined above. It is an immediate consequence of the inequalities (3), (4) and (5) 
that solutions of the form z=q(x, y), defined for all x, y, are necessarily of the 
form z=ax-+by-+c; and more generally that connected complete solution 
surfaces whose normal directions omit a neighborhood of a fixed direction are 
necessarily planes. In the minimal surface case, the former result was first proven 
by BERNSTEIN [5], and the latter was conjectured by NIRENBERG and proven 
by OssERMAN [9]. For non-parametric minimal surfaces, estimates of the second 
derivatives and of the curvature were originally given by HEINz [10], and the 
inequality (4) was given by Hopr [//]. An elementary proof of the Hopr ine- 
quality was given by Nirscue [12]. A generalization of the HEINz estimates to 
a class of non-parametric variational problems has been given by the author [4]. 

Our methods permit the formulation of a problem in partial differential 
equations which is a generalization of the parametric variational problem (1) 
in the same sense that the quasi-linear equation a(y,, Py) Pry +20 (Px; Py) Pry + 
C(Y,, Py) Pyy =O Is a generalization of the Euler equation of the non-parametric 
variational problem 6 f F(p,, y,) dx dy =0. 

If (A(u, v), B(w, v)) is a conservation law of a non-parametric variational 
problem 6 f F(u, v) dx dy =0, then Bdx — Ady is an exact differential on each 
solution surface z=q@(x, y). LOEWNER [6] has shown that the conservation laws 
of a non-parametric variational problem of this type are the solutions of a certain 
first-order linear system. Now let S be an arbitrary solution surface of a regular 
parametric variational problem of the type of (1). In some neighborhood of each 
point g of S, S can be represented in the form z — (x, y) for a suitable coordinate 
system x, y, z; and it follows that the function (x, y) is a solution of a non- 
parametric variational problem of the above type. The system of LOEWNER for 
the conservation laws of this problem is a (non-uniformly) elliptic system and 
hence in the neighborhood of a given point of the w, v plane has infinitely many 
solutions. Thus in a neighborhood of the point g of S there exist many exact 
differentials which arise from conservation laws of the local non-parametric 
problem. We shall investigate the global existence of such exact differentials, 
and we shall show that they arise globally from certain 1-forms on the surface 2. 
These 1-forms are a generalization to the parametric case, of the conservation 
laws of the non-parametric problems. In the case of regular parametric varia- 
tional problems, we shall characterize the generalized conservation laws explicitly 
in terms of solutions of a single uniformly elliptic system. In the minimal surface 
case it will be shown that (a) the conservation laws giving rise to exact differentials 
on minimal surfaces whose spherical images omit a neighborhood are in 1—1 
correspondence with the class of analytic functions regular in a disc; (b) the 
conservation laws giving rise to exact differentials on minimal surfaces whose 
spherical images omit a single point are in 1—1 correspondence with the class 
of entire analytic functions; and (c) the conservation laws giving rise to exact 


13* 
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differentials on minimal surfaces whose spherical images cover the sphere are in 
4—1 correspondence with the set of complex polynomials of degree =2. 


II. Assume that £(p) satisfies conditions (i) and (ii); then as indicated above, | 
the surface Y defined by Y(p) =1 is a closed surface, star-shaped about the 
origin, and the distance R(N) from the origin to the point of ' in the direction | 
of the unit vector N is given by R(N) =1/f(N). Let S be an orientable surface | 
of class C2 in Euclidean 3-space, and assume that one of the normal directions | 
to S has been chosen to be the positive normal direction. Let g be an arbitrary ' 
point of S, and let N(g) be the unit normal vector to S at g. The mapping ; 
p(q) =R(N)N is then a mapping of S into 2 which is a generalization of the » 
spherical mapping into the unit sphere. This mapping will be called the “normal | 
mapping’. Associated with a differentiable mapping of a manifold A into a\ 
manifold B, there is a trivial linear mapping of the differential forms on B into) 
those on A; 1.¢., if %,..+, Xm» Viy «++» Vm are the local coordinates of corresponding { 
points of B and A respectively, then dx; is mapped to 2(dx;/dy,) dy; under the > 
trivial mapping. We now show that the normal mapping of S into Y determines ; 
a non-trivial linear mapping ® of the differential forms on 2 into those on S. 
Assume that S and & lie in different spaces X and Y respectively, but that 
directions in X and Y have been identified so that the normal mapping is well 
defined. In order to define ®, we define specific systems of coordinate patches 
and local coordinates on S and 2. We take for the coordinate patches of 2 the: 
open regions which lie on one side of planes through the origin. For the coordinate: 
patches of S, we take open regions of S whose images under the normal mapping: 
lie in coordinate patches of 2. A coordinate patch V of S whose image under the: 
normal mapping lies in a coordinate patch U of & will be said to “correspond” 
to U. Thus, let o, t, @ be a coordinate system in Y, and let x, y, z be a coordinate* 
system for X, with the x, y, z axes parallel respectively to the o, t, w axes. Let 
U be the coordinate patch of 2’ for which m<0. A coordinate patch V of S: 
which corresponds to U is thus a region of S in which the normals make an acute: 
angle with the negative z-axis. In V, therefore, S can be represented in the form 
z=@(x, y). We take x, y to be the local coordinates in V; and in U we take as 
local coordinates uw, v, where 


(6 weet, gout) 
@ @ 


Equations (6) define a 1—1 mapping of U onto the wu, v plane; since if o, T, ae 
are the coordinates of a point P of U, u,v are the o, t coordinates of the point 
obtained by projecting P through the origin into the plane w =—1. Moreover’ 
the Jacobians 0(t, w)/4(w,v), d(@, o)/A(u, v), O(c, t)/@(u, v) do not simultane- 
ously vanish, since conditions (ji), (ii) imply that the outer normal to Y always 
makes an acute angle with the radial direction; i.e, =o f, 4 tf£,+0f,>0: 


Putting F(u, v) =—1/m and u=g,, v= Py, We see that in V the normal mappin 
has the form 


(7) o=4u/F, t=vF, w—iF. 


Let 6, T,@ be a system of coordinates for Y, obtained from the o, t, w coordinate 


system by a rotation defined by a rotation matrix A = (a;,;). Let U be the coordi: 
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nate patch defined by ®<0, and assume that UXU +0. Noting that the inverse 
of A equals the transpose of A, the relations between yw, v and the local coordi- 
nates u,v in U are easily seen to be 


7 HU th2—%3 FG — Gate. ay5 
(8) ‘ A : A 
ie Ay UTA, ely fea yp UT Ay 20 —A 9 
A A 
where A =— a3, U— 390+ 433, A =— a,,%— ay30-+a33. The region of the u,v 


plane corresponding to UU is the half plane A> 0, since in UNU, w, O<0 and 
@=Aw. A calculation shows also that 4 =1/A, and that @(u, v)/6(a, 3) =A. 
Let V be a coordinate patch of S which corresponds to U and intersects V. Let 
Xx, y, z be the coordinate system for X obtained from the coordinate system %, y, 2 

by the rotation with matrix A. %, 7 are then local coordinates in V, and in V 
S can be Jigs in the form z= A (x, y). The relations between x, y, y(%, y) 


X= 1% + Ay) + 431 (X, ¥) 
(9) = 2% + Ay2V + 4309 (%, ¥) 


P(X, VY) = 43% + Ag3 V + Aga P (%, ¥). 


We now put «=9@,, v =q,, and calculate %, v according to (8). It then follows 
from (9) that #7 =@;, v= ;. A calculation now shows that in VoV, 
A? xz = U5, AP ys = — U5 
10 
( ) A= vy, A? V5 = Uy. 


We now define the ®-mapping for 1-forms. Let / be a 1-form on Soll nas has 
then given by an expression = 8 (u, v) du +A (u,v) dv; and in U, by an expression 
h = 8 (i, 0) dui + A(u,0) dd. In UU, 8, A and #, A must satisfy 

(11) 8=0u,+Arz, A=0uz+Ar;. 


@h is now defined in any coordinate patch V, corresponding to U, by 


(12) Dh =F[A (u(x, y), v(x, y)) dx — B (u(x, y), v(x, y)) dy]. 


This definition is consistent since (8) implies that @/a =A =F/F, and (10) and (11) 
imply that in V0V 

Oh = F2(A dx — Ody) = F? (A xz — Py5) dx + (A 45 — By5) dy] 
03) = F2(Adx — dy). 
Similarly, if f is a quadratic differential form given in U by 
(14) f=Adw+2Bdudv+Cdv?, 
@} is defined in each coordinate patch V of S by 


(15) @j =F4(C dx? — 2Bdxdy + Ady). 
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The consistency of this definition follows as in the case of 1-forms. The ®-mapping 
is obviously linear and may be defined analogously for differential forms of 
higher order. 

If f is a differential form on », we shall denote f considered as a differential 
form on S (i.e., the image of f under the trivial mapping) by the same symbol. 
The following theorem now gives the connection between the ®-mapping and the 
variational problem (1). 

Theorem 1. Let f denote the second fundamental form of X; then a surface S 
of class C2 is a solution of the Euler equations of the variational problem (1) vf and 
only if, at every point g of S, 


(16) {|Pf = — R*(N(q)) K(Q), 
where K(q) is the Gauss curvature of S at q. 


Proof. Let V be a coordinate patch of S corresponding to a coordinate patch 
U of X. Assume that S is given in parametric form with parameters s, ¢ and 
that a region V* of the s,¢ plane corresponds to the coordinate patch V of S. 
Let dA(q) denote the element of area on S at the point g. Then 
(17) Sf Fp) ds dt = [1/R (N(q))] 4A(Q) = eee v)dxdy 

Ve yA ” 

where V’ denotes the projection of V onto the plane of its local coordinates x, y. 
It now follows that S is a solution in V of the Euler equations of (4) if and only 
if the function w(x, y) which represents S in V is a solution of the Euler equation 
of the non-parametric variational problem 6 ff F (u,v) dx dy =0; .e., if and only if 


(1 8) ie uPrx ate ar, v Pry af pe v Pyy = 0. 


Let L(g) denote the expression F,,, 9, ,+2F,y Pry tio Gy, A calculation then 
shows that 


Fi, du? + 2F,,dudv + F,,dv® = L(@) (9,422 + 29,,dxdy + Qyy dy?) + 


(19) 2 
a (p vy Prx Py 9) (F, v a x*® — Pd ips v dxd yy = Die u d yy?) : 


It will be a consequence of Lemma 1, proved below, that in U the second funda- 
mental form f is given by 


(20) f=A(u, v) (F,, av? + 2F,,dudu+F,, dv?) 


vv 


where A(w, v) is a positive function. (Alternatively, equation (20) may be obtained 
by a calculation showing in addition that A = (1/F) (F? + F?+ (F—uF, vF,)?) ya 
Since W/F=R, where W=(1+1?+-v2)), the result now follows from equations 
(15), (19) and (20). In the minimal surface case, (16) reduces to the well known 
formula K(g) =—dso/ds;; since X is then the unit sphere, f =ds? is the Euclidean 
metric on 2 and @f =ds? is the Euclidean metric on S. 

If #(p) is regular, the second fundamental form f is positive definite and 
clearly so is ®/. We may therefore take /, ®f to be Riemannian metrics on Y 


and S respectively. Theorem 1 now implies that the normal mapping of S into X' | 


is conformal in the metrics ®f, f. 


| 
| 
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Lemma 1. Let & be a surface of class C?, represented parametrically in the form 
(u,v) =(o(u, v), t(u, v), wu, v)). Let r*=Pr be a projective transformation of 
the space. x*(u, v) =Pr(u, v) is then a parametric representation of the surface 5* 
which is the image of X under P. The conclusion is that the second fundamental 

forms of X' and &* are proportional at corresponding points. 


Proof. The most general projective transformation of o, t, w space is given by 


o* = (4/6) (Rio +FAet+PR,o0+BfP,), 
| (21) t* = (1/0) (Bio+Ait+h,o0+P,), 
w* = (1/6) (Rio + Bet + Ryo 4+ B,,), 


where 6=P,,6+P,.t+P,0+P,,,and P=(P,;) isa4 x 4 matrix with det (B;) =1. 

Note that we are using the symbol P ambiguously to denote both the projective 

transformation and the 4 x 4 matrix. The precise meaning will be clear from the 

context. Let z, 2* denote the 4-vectors (a, tT, w, 1), (o*, T*, w*, 1) respectively. 

It then follows that z* =(1/6) Pz. Let Ldu®?+2Mdudv4+N dv?, L*dwt+ 

2M*du dv + N*dv? denote the second fundamental forms of 2 X* respectively, 
A 


and let B denote the 4 x 4 matrix whose rows are the 4-vectors A, B, C, D. Then 


Cc 
D 
ayy 
* * Bie 
at, Ox By det eo 
2* 
(1/0) Peis = 2(1/0) x Te ipl (1/6) uu Pz 
— et | (tl8) Pee + 10a Pe 
(1/6) Pz, + (1/0), Pz 
(22) (1/6) Pz 
Pix Cun 
ies z 
= (1/6)tdet |~ | = (4/6)4 det (B,,) det | *" 
(nfo) det |" | = (/8)* det (B,) det | 
IP gs z 


cna (1/0)* (Vip “Ty X ty) S 


Similar equations hold for x2, - vf rf and vj, - ty % ; thus L*/L =M*/M =N*/N= 
(1/6)4 (|v, x v,|/| r* x r*|) =z (wu, v). The step needed in the proof of Theorem 1 
follows since the surface w =F(u, v) is the image under a projective transforma- 
tion of the neighborhood U of 2. This lemma is more elaborate than is necessary 
for the proof of Theorem 1; however, we shall need it again later. 


III. The parametric variational problems considered above may be gener- 
alized to problems in partial differential equations. Let 2, be the unit sphere, 
and let f be a quadratic differential form on 2. Let S be a surface of class C?, 
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and define coordinate patches and local coordinates on S and 2; as before. We 
shall say that S is a solution of the problem associated with / if at each point of S 


(23) [Dj =—K 


where K is the Gauss curvature of S. Thus let U be a coordinate patch of 2), 
and let V be a coordinate patch of S which corresponds to U. It may be readily 
verified in this case that F(w, v) has the same form in every coordinate patch and 


is given by 


(24) F(u,v) = —1jo=(1+w+v0)=W. 
In U, let / be given by 
(25) f =a(u, v) du? + 2b(u,v) dudv- c(u, v) dv?. 


@f is then given in V by 
(26) Of =W!(cdx?— 2bdxdy+ady?’). 


A calculation similar to (19) now shows that the part of S lying in V is a solution 
if and only if the function g(x, y) representing S in V satisfies the quasi-linear 
equation 


(27) ADae +20 Dey + € Pyy = 0. 


If { is a positive definite quadratic form, then @f is positive definite also, and f, 
@f may be taken to be Riemannian metrics on 2, S respectively. Again, if S 
is a solution, the spherical mapping is conformal in the metrics ®/, 7. The follow- 
ing problem suggests itself. Given a Riemannian metric dsj on the unit sphere, 
and a closed Jordan curve in space; does there exist a surface S which spans the 
curve such that ds3/® dst —— K? 


IV. If two Riemannian metrics are defined on the same two-dimensional mani- 
fold, the infinitesimal circle at a given point in one metric is an infinitesimal 
ellipse in the other metric. If the ratio 6 of the major axis to the minor axis is 
uniformly bounded for all points of the manifold, the metrics are said to be 
quasi-conformally related. If *,, x, are local variables in a neighborhood of the 
manifold and the metrics are given by a, dxj+2b,dx,dx,+c,dx3, agdx+ 
2b, dx, d%_+ Cy dx3, then a formal calculation shows that 


(28) (0 | 1/0) E il e Ay Coty C6, — 2D, by 


1 
ghia) 2 (ac, —b2)} (ay cy—03)! 

The quantitiy E is said to be the eccentricity of one metric with respect to the 
other, and clearly the metrics are quasi-conformally related if and only if 
SupE<oo.* If the manifold is compact and the metrics are continuous, then 
clearly the metrics are necessarily quasi-conformally related. In particular any 
continuous metric dsj on the unit sphere is quasi-conformally related to the 
Euclidean metric. In terms of the local coordinates “, v in a coordinate patch U 
of the unit sphere 2}, the Euclidean metric is given by ds? = (1/W4) (4+ v?) duz— 
2uvdudv + (1+?) dv*]. It follows therefore that if f is a continuous, positive 
definite quadratic form on 2; which is given in U by f =a du2+2bdudv+cd UF: 


* The criterion (28) is due to Finn. Fora geometrical interpretation see [7]. 
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IIA 


(29) pakaee a(1+u?)+c(1+v?)+2u0b ise 


2 (ac—b2)* (1+ u2+v2)3 


| This condition on the functions a, b, c is precisely the condition that equation (27) 
be an equation of minimal surface type. 
Now let 2’ be the surface #(p) =1, where Y(p) satisfies conditions (i) and (ii). 
If we map 2 onto the unit sphere X, by projection through the origin, each 
coordinate patch U of 2 is mapped onto a coordinate patch U’ of 2, and corre- 
sponding points of U and U’ have the same local coordinates u,v. We may 
transplant the second fundamental form of 2’ to 2, by means of the trivial map- 
ping, and in terms of the local coordinates wu, v it will be given by the same ex- 
pression whether regarded as a quadratic form on & or as a quadratic form on 3}. 
Thus the variational problem 6 { Y(p) ds dt =0 may be taken to be a special 
case of the problems discussed in Section III. If Y(p) is regular, then the second 
fundamental form of 2 is positive definite and, considered as a metric on 3), 
is quasi-conformally related to the Euclidean metric. Equation (27) coincides in 
this case with equation (18) and, by the above remarks, is an equation of minimal 
surface type. According to the results of Section II, equation (18) is the Euler 
equation of a non-parametric variational problem 6 f F(u, v) dx dy =0, where 
the surface w =F (u, v) is obtained from the neighborhood U of the surface »’ 
by the projective transformation u =—o/w, v=—t/w, F(u, v) =—1/w. Equa- 
tion (18) may therefore be written in the form (F,),+ (F,),=0, and therefore the 
pair (F,, F,) is a conservation law of equation (18). . 

We now show that if Y(p) is regular, then equation (18) satisfies the hypo- 
thesis of the theorems of FINN and of the theorem of the author which are 
mentioned in the Introduction. The theorems of FINN apply to equations of the 
type of equation (27) which are of minimal surface type and have a bounded 
conservation law. The theorem of the author on the other hand applies to the 
Euler equations of non-parametric variational problems of the above type, which 
are of minimal surface type, and which satisfy the condition that the quantity 
|F—uF,—vF,| is bounded. That equation (18) has these additional properties 
may be seen immediately from the following observation. If # is a point of U 
with local coordinates uw, v, then the distance 7(f) from the origin to the tangent 
plane to J at # is given by 


(0) r(p) =1/(F2 + F2 + (F—wF,— vf). 


If Y(p) is regular, then is a closed convex surface containing the origin, and 
y(p) is clearly bounded above and below by positive constants. The result now 
follows from (30). We have now proved the following theorem. 


Theorem 2. Let S be a solution surface of the form z=q@(x, y) of the para- 
metric variational problem 6 f £(p) ds dt=0. Assume that J(p) ts of class C* for 
all p=-0, and that J(p) is positive homogeneous of degree 1, positive definite, and 
vegular. The Euler equation of the non-parametric variational problem 6 fF (u,v) dxdy 
satisfied by (x, y) is then an equation of minimal surface type, having a bounded 
conservation law and satisfying the condition that |F—uF,—vF,| ts bounded. 


190 H. B. JENKINS: 


V. We return now to the study of parametric variational problems satisfying 
conditions (i) and (ii). 2% is a closed star-shaped surface of class C?+o and fis they: 
second fundamental form of X. Let S be a solution surface, and let h be a 1-form 
on X. It is natural to ask if there are conditions on / which will imply that Oh | 
is an exact differential on S. Again let U be a coordinate patch on 2 with local 
variables u, v. Let V be a coordinate patch on S with local variables x, y, which 
corresponds to U. Let / be given in U by 


(31) h=#(u,v)du+A(u, v) dv. 
@h is then given in V by 
(32) Oh =F2(Adx—#dy). 


Thus ®h is exact in V if (F?0),+ (F#.A), =0; 7.¢., if F?0, F2 are the components 
of a conservation law of the quasi-linear equation (18) which is satisfied by the 
function w(x, y) representing S in V. According to the theory of LoOEwNER [6], 
this will be the case if 2%, F2. satisfy the linear system 


| AUS oS) SE Ae Zhe 


(33) 2F,4(F20),, = Ful E20), + (FA), |. 


If this condition holds in all coordinate patches of 2’, then ®/ is an exact differ- 
ential on all of S. In the following lemma we show that with reference to a fixed 
coordinate system in X and Y, the images under the ®-mapping of the differ- 
entials of the coordinate functions on » are exact differentials on S and that the 
differentials of the coordinate functions on S are the images under the ®-mapping 
of differentials on &. 


Lemma 2. Let o,t,@; x, y, z be coordinate systems for Y and X respectively 
with the o, t,@ axes parallel respectively to the x, y,z axes. It then follows that 
(a) dé = @(dt), dn =®(—do), dy =@ (dw) are exact differentials on S and that 
(b) dx = ®(t dw —w dt), dy =®w do —o da), and dz =®(o dt—t do). 

Proof of part (a). Let V be an arbitrary coordinate patch of S, and 
let U be a coordinate patch of X' to which V corresponds. Let 6, T, @ be a co- 
ordinate system for Y such that U is defined by ®©<0, and let %, y, Z be a parallel 
coordinate system for X. X, are then the local coordinates in V. do, at, dw 
are linear combinations of dd, dt,d@; hence it suffices to show that D(do), 
(dt), O(d@) are exact in V. In terms of the local coordinates a, in U, however, 


dé = (1/F)2[F — uF) du — iF] 
(34) dt = (1/F)?[— dR, du + (F —dF) dd] 
do = (1/F)? [Fda + F.dd), 
hence 
O(— dé) = uF, dx + (F —aF,) dy 
(35) P(dt) = (F —0F) dx +0F, dy 


Parametric variational problems 191 


This proves part (a), since it is well known that —iaF,, F—aF.: F—oF, —ivE 


us uU? 


) and F,, F; are the components of conservation laws of the equation 


36) lor eae F392; + Bs G5 =0 


We remark that the exaetaecs of d&, dn, dy can be inferred from the parametric 
form of the Euler equations and hence is well known. Our method has the 
/advantage that it yields quantitative information about the functions &, x, 
which arise by integration. 


Proof of part (b). Let V, U: %, y, 2; 6, T, @; u,v, be defined as in the proof 
of part (a). Note that the assertion of part (b) may be expressed in vector notation 
as dr —@®(RxXdR), where r, R denote the position vectors of points of S 

‘and & respectively. In V it may be readily verified by calculation that 
(dxX%=@O(tdo—Odt), dy=P(O@do—Gdo), d?=—@O(6dt—Tdc). The result 
now follows from the geometrical invariance of dr and Rx dR under rotation. 


In Lemmas 3 and 4, we continue to assume that d&, dy, and dy are obtained 
via the ®-mapping from the differentials of the coordinate functions on SY with 
respect to a fixed coordinate system 6, T,@ for the space Y. We shall assume 
that S is a simply connected solution surface so that the differentials dé, dy, dy 
give rise, by integration, to single-valued functions &, 7, y on S. We now derive 
some properties of the mapping of S into the plane which is defined by the func- 
tions &, 7. 


Lemma 3. Let g be a point of S, and let p(q) denote the image of gq on X under 
the normal mapping ; then (a) €,n 1s L—I1 in a neighborhood of q provided that the 
w-component of the normal vector to X at p(q) does not vanish; and (b) 


(37) cos? #/R? < (dE/ds)? + (dylds)? < 1/R* sin? gp, 


where ds is the Euclidean element of length on S, 3 ts the angle between the tangent 
plane to X at p(q) and the o, t plane, and is the angle between the tangent plane to 
»' at p(q) and the radial direction. 


Proof of part (a). Let V be a coordinate patch of S containing g, and such 
that the plane of the local coordinates x,y is parallel to the tangent plane to 
S at g. Let U be a coordinate patch of 2, with local coordinates #, 0, to which 
V corresponds. It follows that 0(&, 7)/é(%, ¥) =F*[8(o, t)/@ (a, 3)] 0; since 
6(o, t)/0(%,v) is the w-component of R;xR;z. R again denotes the position 
vector of points on 2. 


Proof of part (b). If /, g are quadratic forms on &, it is clear from the defini- 
tion that in corresponding coordinate patches V, U of S and &, respectively, the 
ratio @f/®¢g at a given point of V in the direction dy/dx =m has the same value 
as the ratio f/g at the normal image point on U in the direction dv/du =— 1/m. 
The result now follows since d&+d7?= @O(do?+dr?), ds?=@O(R x dR)?, and 
clearly 
(38) cos? #/R? < (do® + dt?) /(RX dR)? < 1/R* sin? ep. 
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In a neighborhood of a point p of X at which the normal vector has a non- 
vanishing w-component, the w-coordinate on » may be taken to be a function 
of o, t. According to the preceding lemma, in a neighborhood of a point q, whose 
normal image on » is ~, &, 4 may be taken as independent variables. Parts (a) 
and (b) of the following lemma are due to Rapo. 

Lemma 4. In a neighborhood of a point q, where the &,n mapping ts I—I; 
(a) z =o, %, =7, (b) 2(&, 7) 7s a solution of the Euler equation of the non-parametric 


» &n 
variational problem 6 { w(o,t) d& dn =0; and (c) K(q) = we O(a, t)/0(E, ), where 
K(q) is the Gauss curvature of S at q. 

Proof of part (a). Observe that in any coordinate patch V of S, which 
corresponds to a coordinate patch U of X; dy =0 in the direction dy/d% =o;/0;, 
and that d& =O in the direction dy/dx =t;/t;. The result now follows from 
Lemma 2. 

Proof of part (b). To prove part (b), note that the Euler equation of the 
variational problem involved may be written in the form (w,):+(@,),=0. The 
result now follows from Lemma 2; since in any coordinate patch V of S, corre- 
sponding to a coordinate patch Taf 2, 


(39) w, dé —w,dn = F®(w;d% — wz dy) = dy. 


Proof of part (c). Let V be a coordinate patch of S containing g and 
corresponding to a coordinate patch U of 2, then 


0) KM =(5) aay Ce) Leen aaa] een 
but £ 

\4f a(é, T fe NA 
wt (sa) lace ay/-aeeay] = Ge) = ae 


In the case of non-parametric variational problems of the type 6/F (u,v) dxdx=0 
the addition of a linear function of w, v to the integrand F(z, v) leaves the Euler 
equation and hence also the class of solutions unchanged. An analogous situation 
exists in the case of the parametric variational problems under consideration. 
Since the integrand ¥(p) of the variational problem (1) is completely defined 
by its values on the unit vectors, £(p) is uniquely determined by the surface 2. 
We now show that 2’ is not uniquely associated with the class of solution surfaces 
but may be chosen according to the following lemma. 


Lemma 5. Let P be a projective transformation of the space Y which maps 
each line through the origin onto itself, preserving orientation. Assume also that P 
maps every point of & to a finite point. Let &* denote the images of S under the 
transformation P; i.e., 2* =P. The conclusion is that the class of solutions is 
unchanged if X' is replaced with *, * 


Proof. Let o, t,@ be an arbitrary coordinate system for Y, and let U, U* 
be the coordinate patches of 2’, XY* respectively for which w<0. The hypotheses 
on P imply that PU=U* and that points of U, U* which correspond under P 


x it can be readily verified that replacing Y by S* has the effect of adding linear 
functions to the integrands of the local non-parametric variational problems. 
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have the same local coordinates u,v. Let f, f* denote the second fundamental 
forms of 2’, X* respectively; then in U according to Lemma 1, f=y(u, v) f*. 
| Let S be a solution surface of the problem associated with 2; then by Theorem 1, 
| {/®f=— R‘K. Let V be a coordinate patch of S corresponding to U, and hence 
also to U*, and let F*, R* be the quantities associated with Y* which are analo- 
_ gous to F, R (1.e., F* (u,v) is the negative reciprocal of the w-coordinate of the 
point of U* whose local coordinates are w, v) pihensinceyk Li == RAR, 


(42) eT (EGP) (HO) = — ReEK 


The hypotheses of Lemma 5 can be weakened in the following sense. Let 4” be 
a subset of 2, and let P be a projective transformation which satisfies the hypo- 
theses of Lemma 5, except that instead of requiring that P map all points of X 
to finite points, we require only that P maps all points of X —.¥ to finite points. 
The class of solution surfaces whose images under the normal mapping omits the 
set VW of 2’ is then unchanged if » is replaced with 2* = PS. 


VI. In Sections II and IV we have assumed that £(p) is positive homogene- 
ous and positive definite. We now assume in addition that £(p) is regular; 7.e., 
»' is a closed surface of positive curvature. 


Lemma 6. Let x, y, z be a fixed coordinate system for X, and let F denote the 
class of solution surfaces of the variational problem associated with SX whose normal 
directions ont the positive z direction. Assume that X' is a surface of positive curva- 
ture and hence 1s convex. Leto, t, w be a coordinate system for Y, with axes parallel 
respectively to the x,y,z axes. Let py be the point of intersection of &X' with the 
positive w-axis. There then exists a projective transformation P with the following 
properties: (a) P maps po to a point at infinity, P maps all points of X'— po to 
finite points, and otherwise P satisfies the hypotheses of Lemma 5; (b) the class 
1s unchanged if X' 1s replaced with XY* = PD’; (c) the surface Y* = PL 1s represent- 
able in the form w =G(o, t) where G ts defined for all o, t; and (a) the Euler equa- 
tion of the variational problem 6 f G(o,t)d&dy=0 is uniformly elliptic in the 
entire o, t plane. 


Corollary. Let S be a simply connected solution surface whose normal directions 
at each point make an angle =P with the positive z direction. Let do, dt be the 
differentials of the coordinate functions on X*. The functions &, yn arising by integra- 
tion from the exact differentials dé =@(dt), dn =®(—do) then define a mapping 
of S into the plane which is, (a) locally 1—I1 everywhere on S; and (b) has the 
property that (dé/ds)?+ (dn/ds)?=c? at every point of S, where ds is the Euclidean 
element of length on S, and c ts a constant depending only on B and &. 


Proof of the corollary. The normal image of S lies in a compact subset of 2*. 
»* however is a surface of positive curvature of the form w =G(g, t), where G 
is defined for all values of c, t. The inclination of the tangent plane to 2™* is 
therefore bounded on compact subsets. The corollary is therefore an immediate 
consequence of Lemma 5, and parts (a) and (b) of Lemma 3. 


Proof of Lemma 6. Let 7 denote the tangent plane to 2 at 9, and let H 
denote the plane through the origin parallel to T. Let 6, T, © be a coordinate 
system for Y such that the 6, T plane coincides with H, and T intersects the 


194 H, B. JENKINS: 


positive @-axis. We now define the projective transformation P in terms of the 
G, T, ® coordinates as follows: 


- ee ee ee ee 
1—a/h’ 1—@/h —Oa/h 


(43) Cries 


where h is the perpendicular distance from the origin to the plane T. From (43) 
it follows immediately that P maps each line through the origin onto itself 
preserving orientation; that P maps fy to infinity, and all points of '— pp to finite 
points; and that P maps the plane 7 to infinity and leaves the plane H pointwise 
fixed. The equations (43) 
also imply that the image 
under P of a line L inter- 
secting 7 in the point Q’ 
and H in the point Q”, 
is the line L’ through Q” 
and parallel to the line 
through the origin and Q”. . 
The surface 2* = P22’ and 
also the function G(o, T) 
can now be constructed 
in the fashion of projec- 
tive geometry (see Fig. 1). 
Let Q be an arbitrary 
point in the o, Tt plane, — 
and let L, be the line 
through Q, parallel to the 
j-axis. Let Q, be the point 
where JL, intersects the 
plane H, and let L, denote 
the line through Q, and fp. 
According to the above remarks, then, the line L, is mapped onto the line L, 
under the transformation P. Let p be the point where the segment fy Q, intersects 
2’, and let L, be the line through the origin and #. Since L, is mapped onto L, 
by P, and L, is mapped onto itself, it follows that the image of # under P is the 
point * in which J, intersects L,. If o, t are the coordinates of Q, then G(g, 7) 
is the height of £* above the a, t plane. This proves parts (a) and (c) of Lemma 6, 
and part (b) follows directly from Lemma 5. 


|@ 


Fig. 1. Construction of 2* 


Proof of part (d). We first observe that the ellipticity of the Euler equation 
of the variational problem 6 f G(o, t) dé dy =0 is equivalent to the condition 
that 2* is a surface of positive curvature and hence follows from Lemma 1. 
The condition of uniform ellipticity is equivalent to the condition that the second 
fundamental form of X*, considered as a Riemannian metric over the o, t plane, 
is quasi-conformally related to the metric do? dt?. According to equation (28) 
this will be the case if and only if there is a constant & such that 
(44) __ Geet Gow 


— <> ROO} 
(Goo Gz,—G? 1)? 
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We note next that if R denotes the position vector of a surface and if the surface 
has the property that the angle between the outer normal vector and R is acute 
at each point, then (R x dR)? is a positive definite quadratic form and may be 
taken as a Riemannian metric on the surface. Furthermore, if R* is the position 
vector of another such surface, then the mapping of one surface on the other by 
projection through the origin is conformal in the metrics (R x dR), (R* x dR*)?. 
This follows since such a mapping has the form R* =/R where A is a scalar 
function, and hence R* xdR*=ARx(AdR + (dd) R) =2 (RX dR). Now let 
R, R* denote the position 
vectors of 2’, ** respective- 
ly. We shall show that the 
metric (R* xd R*)? is quasi- 
conformally related to the 
metric do*+ dt. Part (d) 
of Lemma 6 then follows 
from the transitivity of the 
relation of quasi-conform- 
ality. To see this, observe 
that by the above re- 
marks, the transformation 
P maps 2* conformally 
onto 2 —, with respect 
to the metrics (R* x dR*)?, 
(R x dR)?. The metric 
(Rx dR)2 on 2 is neces- 
sarily quasi-conformally re- 
lated to the second funda- 
mental form of 2’considered 
as a metric, by the com- 
pactness of 2. On the other 
hand, Lemma1 imphes that Fig. 2 
the transformation P is con- 
formal with respect to the second fundamental forms considered as metrics on 
+, »*. It remains to show that the metric (R* x dR*)? is quasi-conformally 
related to the metric do?+dt?. R*¥= (o, t, G(o, t)), and it now follows from a 
calculation using equation (28) that what must be shown is that 


Ww 


R*?+ (G—oG,—tG,)*+ (6G,—tGe)?* 
(45) RS (G=6G,= 76a 


where R* = |R*|. Let Q=(o?4+ 72)!, and put Z =|G—oG,—1G,| =|G—NGo|. 
Let L,, L,, L, be defined as in Fig. 1, and let y denote the plane containing the 
lines L,, L,, Lz. Let L4 be the tangent line to 2’ at #, in the plane y (see Fig. 2). 
Let Q,, 0, be the points where L, intersects the plane H and T respectively. Let 
L, denote the line through Q, and #*; then under the transformation P, L, 1s 
mapped onto L;. L; is therefore the tangent line to 2* at p* in the plane y. 
Since y contains the w-axis, the slope of L, is the quantitiy Gg. It follows that 
the quantitiy Z=|G—QGo| is equal to the distance from the origin to the 


cate <A, 
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point where L, intersects the negative w-axis. Let o% be the angle between L,, 
and the w-axis, and let « be the angle between L, and L;. Application of thelaw 
of sines shows that Z/R* =sina,/sing,. The line L; is parallel to the line Ly 
through the origin and Q,; therefore « is equal to the angle between L; and Lg. 
Let I’ be the curve of intersection of X’ and the plane y. It may be shown by; 
elementary means that the limit as p approaches fy along J’ of sina/sina, is: 
equal to one less than the order of contact of J’ with its tangent line at fy. Sincet 
the curvature of J” does not vanish, the order of contact is 2, hence the limit in: 
question is 1. From the smoothness properties of it may now be inferred that: 
as a function of a point p on 2, the quantity Z/R* is continuous and positive. : 
Thus there exist positive constants k,, k, such that k}<Z/R*<k,. The left-hand: 
side of (45) may be written in the form R*/Z+Z/R* + (oG,—1G,)?/(R*Z), and: 
the boundedness of the first two terms has now been established. It is sufficient: 
to show that the last term is bounded in the exterior of a compact region of thes 
o,t plane. Let J; denote the level curve of 2* at height /. Since 2* containss 
the origin, for all h=O, Ij is a closed convex curve enclosing the w-axis. The¢ 
curve J; is in fact the projection through the origin of a certain plane section of 2° 
onto the horizontal plane at height #. It may be verified easily that the level: 
curve at infinite height is the projection through the origin of the Dupin indicatrixs 
of X at p, onto the horizontal plane at infinity and hence is an ellipse of finites 
eccentricity with center on the w-axis. It now follows that the angle 6 between: 
VG and the radial direction in the o, t plane is bounded from $2 in the exterior 
of Ij; i.e., there is an angle 6, such that 0<d<6)< 4a in the exterior of J}. 
In the exterior of Jj, therefore, |¢G,—tG,| = 2|VG| sin 6 =Q2Gp tan 6<. 
QGe tan dy, and hence 
(46) (oG,—tG,)? (2Ge tan do)? < (R*+Z)? tan? dy 
eZ, Vie VA aa riety 

Lemma 7. Let S be a simply connected surface of class C?. Let q be a point of S,. 
and let d be the distance along S from q to the boundary of S. Let €, n be differentiable 
functions on S. Assume that (i) the mapping defined by the functions &, yn of S inte 
the plane is locally 1—1 at all points of S and that (ii) (d&/ds)?+ (dy/ds)?=c2s 
where ds 1s the Euclidean element of length on S and c is a positive constant. They 
conclusion 1s that the functions E,n map a subregion of the geodesic circle of radius a: 
about q on S homeomorphically onto the circle (E—& (q))?+ (y—n(g))2< (cd)®. If St 
zs complete, it follows that the functions €,y map S homeomorphically onto thes 
entire E,» plane. 


Proof. Let 7 be the radius of the largest circle in the &,7 plane about the: 
point (&(g), 7 (q)) in which the inverse mapping is 1—1. Suppose that r<cd. 
Let / denote the pre-image of any given radius of the circle (gé—& (q))? + 
(7-7 (q))?<7?. The length of 7 is then necessarily <d since 


(47) x a ((dE/ds)? + (dy|ds)?)' ds =c-length of J. 


Since the mapping is locally 1—1, the inverse mapping is therefore 1—1 in a 
neighborhood of every boundary point of the circle (§—€(¢))?+ (n—n@))2< 
and hence is 1— 1 in some circle of radius >7. A different proof, which is adapta- 
ble to the present case, is given in [4] for the case in which S is a plane region. 
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Let ds?=ado?+2bdodt+cdt* be a Riemannian metric defined over a 
-Yegion D of the o,t plane. Let g=o+i7t, (0) =f(o, t) +7q(c,7); then the 
-mapping 7(@) is a conformal mapping of the Riemannian manifold defined over 
_D by the metric ds? provided #, q is a solution of the Beltrami system with 
coefficients a, b,c. This Beltrami system is given in complex form by 


(48) t= (0) ho 

_ where yw (0) =(a—c+21b)/(a+c+2(ac— b?)#). The Beltrami equation (48) is uni- 
formly elliptic if there is a constant k such that (a+c)/(ac—b?)i<k<oo or 
equivalently if there is a constant jy such that 0< |u|<puyp<1. 


Lemma 8. Let 7(0) be a solution of the Beltrami system (48) in the disc |o| <1. 
Assume that |y|<M, |u|Spo<1 and that (0) satisfies the uniform Hélder 
condition | (0) — 4 (02| SH | e,—05|*. Let o be an interior point of the disc |o| <1. 
Let d, be the distance from @ to the boundary. Then there is a constant k(a, H, |), 
depending only on the quantities indicated, such that d,| 4, (0)|<k(«, H, 9) M. This 
lemma is a special case of a theorem of DouGc.is & NIRENBERG [13] on Schauder 
estimates for elliptic systems. A proof for Beltrami systems is given in [4]. 


At a point where the Jacobian does not vanish, a differentiable mapping 
carries infinitesimal circles into infinitesimal ellipses. If 6 is the ratio of the minor 
axis to the major axis at a given point, then the eccentricity E of the mapping 
at the point is defined to be E=}$(6+1/6). A mapping which is continuously 
differentiable in a domain D is defined to be quasi-conformal in D if its eccen- 
tricity is uniformly bounded there. For a mapping u(x, yy), v(x, y), this definition 
is equivalent to the condition that the Jacobian 0(u, v)/0(x,y) vanishes only at 
isolated points, and that the metric du?-+ dv? is quasi-conformally related to the 
metric dx?-+-dy?. It can be shown that the eccentricity of a mapping given by a 
solution of the Beltrami system 7, = (0), is given by E=(1+|m|?)/(1—|u|?). 
Thus such a mapping is quasi-conformal if the Beltrami system is uniformly 
elliptic. If m(x,¥) is a solution of an elliptic equation agy,,+2b@,,+cC@y,=9, 
ac— b?>0, it can be shown that the eccentricity E of the mapping u(x, y), v(x, y), 
where u=q,, Uv =Q,, Satisfies an inequality ES }((a+c)?/(ac— 6?) —1). Such 
a mapping is therefore quasi-conformal if the equation satisfied by (x,y) is 
uniformly elliptic. 

Lemma 9. Let w=u-+dtv be continuously differentiable with eccentricity SE, 
in a domain A of; the z-plane. Assume that |w|<1 then in any compact subregion 
B of A, w(z) satisfies a uniform Holder condition 


(49) | @ (a) — @ (%)| SH 


where 6 =E,—(Eo—1)!, H is an absolute constant, and d is the distance from B 
to the boundary of A. 

Proofs of the Holder continuity of quasi-conformal mappings have been given 
by many authors. For further discussion, and a proof of Lemma 9 with H=ze, 
see [14]. 

The following lemma is due to WARSCHAWSKI [195]. 

Lemma 10. Let R be a region in the z-plane bounded by a continuously dtffer- 
entiable closed Jordan curve C. Assume that as a function of arc length on C, the 
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angle 0, between the tangent line to C and the x-direction, is Holder continuous with 
constant k and exponent v. Let D be the diameter of C. Let % be an interior pownt — 
of R, and let 0 be the radius of the largest circle about 2 lying entirely im Rk. Lé : 
d—GLB (ra), where x is the distance between any two points of C, and o ts the length ; 
of the shorter arc of C connecting them. Let w=f(z) be the Riemann mapping 
function which maps R conformally and homeomorphically onto the disc |w|<4, 
with f(z) =0. The conclusion is that there are constants py, ly Such that 0<j4S 

| 7’ (2)| SMe, for all z in R, where 1, 4, depend only on k, v, Dip; d: 


We now establish the inequalities mentioned in the Introduction. 


Theorem 3. Let S be a simply connected solution surface of the variational 
problem 6 f J(p)dsdt=0. Assume that J (p)eC?™ for all p+0, and that J (p) 
is positive definite, positive homogeneous of degree 1, and regular. Assume that the 
normal vector at each point of S makes an angle = with the positive z-axis. Let q 
be a point of S, and d be the distance along S from q to the boundary. Let « be the 
angle between the normal at q and the positive z-direction. Then 

=e 2 
(50) igi pyre 
where K(q) ts the Gaussian curvature of S at g, and C(B) 1s a constant which depends 
only on B and the integrand J(p). If S is a surface of the form z= (x,y), then 


2 C 
(51) |K (0, ¥o)| S Gaye 
and ; 
CLW 
(52) % +28 +% Ss, 


where K(X, Vo) 1s the Gaussian curvature of S at x, Vo, v is the radius of the largest 
circle about x, Vo in which p(x,¥) is defined, Wy = (1+ 9% (xo, Yo) +95 (%o» Vo), 
ANd 7, Sq, ty are the derivatives Pyy (Xp, Vo)» Pry (Xo» Vo)» Pyy(Xo» Vo) respectively. 


Corollary. Under the same hypotheses, if S is complete, then S is a plane. 


Proof. Let 2* be the surface constructed according to Lemma 6. The right 
circular cone having vertex at the origin and generators making an angle f with 
the positive w-axis intersects 2* in a closed curve J(3; which bounds a compact 
region D*(f) on X*. By assumption, the image of S on Y* under the normal 
mapping lies in D*(B). Let D(B) be the projection of D*(8) onto the o, t plane. 
Our smoothness assumptions on XY imply that the function G(o, t) is of class 
C?+, uniformly so on compact subsets. Thus it follows that the boundary J"(B) 
of D(B) is a closed convex curve of class C?**. Let &,7 be the functions on S$ 
arising from the differentials of the coordinate functions on ©* by means of the 
@-mapping and Lemmas 2, 5 and 6. It may be assumed that €(¢) =7(¢) =0 
By Lemmas 6 and 7, the functions £, 7 map a subregion of the geodesic circle on 
S of radius d about q¢ homeomorphically onto the circle &24 72< c2d?, where c is 
the constant of the corollary to Lemma 6. According to Lemmas 4 and 6, in 
&+7?<cd®, z=, z,=7, and 2(&, 7) is a solution of the Euler equation of the 
non-parametric variational problem 6 f G(,t)d&dy=0. Thus 2(&,) satisfies 
the equation 


(53) Goo keg + 2Go_ 2 q+ Gre nn =O. 
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By part (d) of Lemma 6, this equation is uniformly elliptic, and according to the 
remarks made above, the values assumed by the derivatives 2 =0,.2, 
the region D (f). 


The corollary follows, since if S is a complete surface, by Lemma 7 the , " 
mapping is a 1—1 mapping of S onto the plane. 0(&) =o(€, 4) —it(, 7) is 
then bounded and quasi-formal in the entire plane. Since Liouville’s theorem 
is valid for quasi-conformal mappings, it follows that o(f) is a constant. The 
normal image of S on 2* is therefore a single point, hence S is a plane. 

By part (c) of Lemma 4 and Lemma 6, the Gaussian curvature K(q) of S at q 
is given by 


(54) K(q) = (R*4) Se) 


, =T lie in 


£=n=0- 


Let H(G) denote the Riemannian manifold defined over the o, t plane by the 
metric dsj =G,,do?+2G,,dodt+G,,dv. Let K(G, z(€,7)) denote the Rie- 
mannian manifold defined over the disc ?+-7?< c?d? by the metric ds3 =G,,d&’— 
2G,,4— dn +G,,dn?. A calculation based on equation (53) shows that 


(55) G,,do?+2G,,dodt+G,,40?=(2,, —Z¢¢%pn) (Gz1d€2—2G,,dé dn +G,,dn?). 


It is well known that z:—%.2,,, vanishes only at isolated points if z(&, 7) is a 
solution of equation (53); hence the function @(¢) =a(&,7)—17(€é, 7) maps 
K (G, z(&, 4)) conformally into H(G), where we have put €=&+7%n. 

In the minimal surface case, 2’ is the unit sphere, and 2* is the paraboloid 
of revolution given by w =}(o?+1?—1). Equation (53) is then the Laplace 
equation, and the projective mapping P of 2X into 2* followed by the parallel 
projection of * onto the plane is just the stereographic projection of 2’ onto the 
o,t plane. The region D(f) is then the circle of radius cot $f about the origin, 
and 9(¢) is a complex analytic function with values in D(f). O(a, t)/0(€, 4) can 
then be estimated by means of Schwarz’s lemma, and the estimate of the curva- 
ture follows from equation (54). This procedure closely parallels OSSERMAN’s 
proof for minimal surfaces. 


In the general case, in order to estimate 0(¢, t)/0(§, 7), we introduce uni- 
formizing variables y(@), y(¢) for the manifolds H(G), K(G, z(&, 7)) respectively. 
The coefficients of the metric ds, are Hélder-continuous; hence by the general 
uniformization theorem there is a conformal homeomorphism ¢(0), with Hélder- 
continuous first derivatives, non-vanishing Jacobian, and such that ¢(0) =0, 
which maps H(G) onto either the disc |t|<1, or the entire ¢-plane. It must be 
the latter, however, since the uniform ellipticity of equation (53) implies that the 
metric ds, is quasi-conformally related to the metric do?+dt?, and the relation 
of quasi-conformality preserves conformal type. Under the mapping ¢(g), the 
region D(f) is mapped onto a region D’(f) in the ¢-plane. The smoothness pro- 
perties of the curve ["(8) and the mapping ¢(g) imply that the boundary Vai (6) 
of D’ (B) is a closed curve of class C**, Let y(t) be the Riemann mapping function 
which maps D’ (f) conformally onto the unit circle | y|<1, with 7 (0) =0. Accord- 
ing to Lemma 10, then, there exist two positive constants f, “, which depend 
only on the Hélder continuity and the geometry of the boundary curve 1” (), 


14* 
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and hence only on the angle f and £(f), such that for all ¢¢ D’ (B), 
(56) (4 S| (dz /4t)| Sue. 


z(o) is a conformal mapping of H(G); hence z(0) satisfies the Beltrami system 
43 = (0) %, Where 
(57) f(Q) = 


1S NG@eg= Ge, F2IG en 
(Gog Cee ais 2IGsa Ger a G2,)?) 


Since equation (53) is uniformly elliptic, 7(@) is quasi-conformal with eccentricity 
<E,)=Sup ((Goo+ Ci Ne (Gy5Gr,—G51)4)), and | (0) | SMy= (Ey—1)3/(Eo+ 1)?<1. 
Since @(t, ¢)/@(0, @) does not vanish in the entire o, t plane, it is bounded above 
and below by positive constants on the compact set D(f). 


Since | 
oo ateca) Vel (4 — le (@) |”) a(t,f) (9.8) ’ 
dy (2 | 
(59) oe <2 (Lael? + laal) =21zel? (1+ 1H). 


where ds is the Euclidean element of length in the a, t plane, it follows that there 
are positive constants m,, mM, N,, N, Such that at all points of D (8), 
dy, 
(60) Mm, =|y¥,| = Me, m S|" smy. 
A conformal mapping y(¢) of A(G, z(€, 7)) must satisfy the Beltrami system 
yz=—u(e())y~. As previously noted, the second derivatives of G(c, t), and 
hence the Beltrami coefficient w(e) are in class C*, and uniformly so on compact 


sets. Thus let H, be the Hélder constant of w(o) on D(8). By the remarks preced- 
ing Lemma 9, @(€) is quasi-conformal with eccentricity <F,, where 


Gre Gaa)® 2 
(61) E, = Sup(576 ee 1) = 2E5 — 4. 


Let a Be aia as then |@(¢)| SM. The Hélder continuity of the map- 
ping o(6) can then be estimated in the disc |¢|<3cd by means of Lemma 9. 
Combining the estimate thus obtained with the Hélder continuity of ju (@) 
see that in the disc |f|<}cd, w(o(¢)) satisfies a Hélder condition 

een a 0; 


(62) | uu (0 (01)) — mw (0 (9) | S A, (HM)* ina 


, 


, we | 


where 6,=E,—(E£j—1)!, and H is the absolute constant of Lemma 9. By the 
general uniformization theorem, there is a conformal homeomorphism y (C) of the - 


submanifold of K(G, 2(€, 7)) defined by |f|<$ed, onto the disc |y|<#ced such 
that y(0) =0. We now show that 


(63) lv; (0)| Sy, 


where K, is a constant depending only on the angle Band Y(p). To see this, put 


y*=(2y)/(cd), C*=(2C)/(cd), and u*(C*) =—w(o(gcde*)). In |¢*)<1;°thent 
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|y*|<1, w*(C*) is uniformly Hélder continuous with constant (HM)*H, and 
exponent «0, and y*(¢*) satisfies the Beltrami system yz, =* (0*) Yee. Since 
ye (C) =y% (C*), the result now follows from Lemma 8. 

It follows from the definitions of y(@) and y(¢) that 7 is a complex analytic 
function of y in the disc |y|<}cd. y,(y) can therefore be estimated at the origin 
by applying Schwarz’s lemma to the function (y— 7 )/(1—7Zpo), where 0) = 
o(0, 0) —77(0, 0), and 7) =7 (09). The result is 


a eee 

Now 

6 8(6,t) _ 2(0,0) a(x, %) aly.) 
a (5, ) 0(%,%) Oy.) a(e,) 


Since 0(y,7)/0(C, 6) =|y-|?(1—|w (0 (C)) |?), it follows from (58), (60), (63) and (64) 
that 

0(4, T) 
66 
a | e(§, 7) 
Now, let y, denote the plane through the w-axis and the point (o(0, 0), 7(0, 0)). 
Let L denote that segment of the line in which y, intersects the o, t plane which 
lies between the point (¢(0, 0), 7(0, 0)) and the boundary of D(f). Let by denote 
the length of L. Under the mapping x(0), the segment L is mapped onto an arc 
connecting the point y) with the boundary of the disc |y|<1. Thus, letting s 
denote arc length in the o, t plane, and using (60), we see that 


_ [x OP lre(0)? Zi ae 
E=n=0 | %o (@0) |? — \ mye a? , 


a. 
(67) t— [nol Sf |S4]ds <n by. 
Ib 


Let L* denote the arc of the curve in which y, intersects 2* which lies over L. 
Let m be the angle between the line from the origin to a point of L* and the 
positive w-axis. It follows by calculation that |ds/dy| =R*(R*/(QGo—G)). In 
the proof of part (d) of Lemma 6, it was established that k, < (QGp — G)/R* Shy, 
where k,, ky are positive constants depending only on the geometry of 2. There- 
fore 


. ds 
(68) by = [|= |dy < (RBA) @ —B), 
B 
where R* (8) = max O*. Combining (54), (66), (67) and (68), and putting 
prs p* * 
Rj = min Of*, we have 
p*es™ 
4 R* (B) 4K, m, |? (x—B)? 
(69) Oe erecta Gems! 


For the case in which S is a surface of the form z=(x,¥), the above results 
are valid with B =42. D(4z) is the region of the o, t plane bounded by the 
convex curve J’(}z) in which 2* cuts the o, t plane, and D* (5) is the region of 2* 
for which G(o, t)<0. For all points (a, t) € D ($2), let b(¢, t) denote the distance 
along the ray through the origin from (¢, t) to the boundary of D(¥x). Let 
by =b(a(0, 0), (0, 0)) and Gy=G(o(0, 0), (0, 0)). Let B= max (—2/G). 


(9,7) €D(3 7) 
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Clearly B is a finite positive number, since (—/G) is a positive continuous func- | 
tion in the interior of D(4z) and takes the values of 1/Gg on the boundary. Thus 


(70) bb = BGy. 


The image under the normal mapping of a point (x,¥v, p(x,y)) on S is the point 
(c, t, G(o, t)) on X*; hence for all (a, t) € D (32) 


(71) W= (1+ 92+ 95)? = (1+ (0/6)? + (1/G)?)* = — (R*/G). 


Now, let (%, Yo, 7 (Xo, Yo)) be the coordinates of the point g on S. Let w= 
(1 +9: (%o, Vo) + 9% (Xo, yo))?, and let v be the radius of the largest circle about 
(Xo, Vo) in which (x,y) is defined. Clearly r<d; hence, combining (54), (66), (67), 
(70) and (71), we have 

a 4K,n, BR*($2) ]2 1 
(72) K(X, Vo) = : mi cRé = We v2 


In order to prove the inequality (52), we shall need the following well known 
lemma. 


Lemma 11. Let S be the surface defined by z=q(x,v) where y(x,y) ts a 
solution of an equation of minimal surface type; 1t.e., an equation of the form 
a(U,v) Py» +26 (u,v) p,, +c (u,v) p,,=0, where the coefficients satisfy the inequality 
(29). The conclusion is that (a) the spherical mapping is quasi-conformal with 
eccentricity <e,, where e, 1s a constant which depends only on the equation; and 
(b) the ratio of the magnitudes of the principal curvatures of S is bounded above and 
below by the positive constants (6,)*, (6,)~* respectively, where 6,=e,+ (ej —1)*. 


Proof. By means of the definitions of Section II, take u,v to be the local | 
coordinates on the hemisphere U which contains the spherical image of S. The — 
spherical mapping is then given by wu =q,, v=q,. The condition of minimal — 
surface type implies that the metric ds?=adu?+2bdudv+cdv? is quasi-— 
conformally related to the Euclidean metric on U, with eccentricity bounded by — 
a constant ey). This same condition also implies that taken as a metric on S, 
P(ds?) =W4(c dx®— 2b dx dy +a dy?) is quasi-conformally related to the Eucli- 
dean metric on S, with eccentricity bounded by the same constant e,. As in 
Sections II, III the spherical mapping is conformal in the metrics @(ds?), ds®. 
The spherical mapping may therefore be obtained by composition from a con- 
formal mapping and two mappings of eccentricity <e,. The eccentricity of such 
a mapping is readily shown to be bounded by e, =2e3 —1. Part (b) of Lemma 11 
is valid for any surface whose spherical mapping is quasi-conformal with eccen- 
tricity <e,. A proof follows easily from the classical formula of differential 
geometry 
(73) JD = OT LT le des) 


where /IJ, IJ, I are the first, second and third fundamental forms, and H, kK 
are the mean curvature and the Gaussian curvature respectively. The normal 
curvature & in a direction making angle g with a principal direction on S is | 
given by 


(74) k =II/I =k, sin? y + ky cos? @, 
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where k, and k, are the principal curvatures. Using (73) and putting 6= 
(max JII/I)/(minIII/I), we see that the eccentricity E of the spherical mapping 
is given by 


(75) E=2 (6+ 1/6) = $ ((Ay/Ra)? + (Ro/Py)?) , 
and the result follows by requiring that E<e,. 


| Now let S be a solution surface of the parametric variational problem under 

consideration of the form z=q(x,y). At a point of S, let , be the principal 

curvature of greatest magnitude. From part (b) of Lemma 11 it follows that 
| %1|S(6,)* |4,|. Thus if k is the normal curvature in any direction, then |k|< 

(Ri)? S (6, |K])#. Now let ky), Ry) denote the normal curvature in the directions 
dy =0, dx =0 respectively; then 


& Pu x = Ty ya 
(76) Ry) = W(1+q?) ’ Ry) W (1+ 2) z 


and it follows that 


(1/W4) [(Pex + Pyy)® + 2 (Psy — Pax Pyy)] SW? (| Rey] + | Agy|)? +2 
< (4W? 6, + 2) |K]. 


K 


(77) 


The inequality (52) now follows by estimating K(%9, Vg) by means of (51). 


VII. In this section we consider the problem of determining all 1-forms on 2’ 
whose ®-images are exact differentials on solution surfaces. We have noted in 
Section V that a 1-form on 2 will have this property if and only if in each co- 
ordinate patch U its coefficients with respect to the local coordinates w, v satisfy 
the system (33). As previously mentioned, such a 1-form is a generalization of a 
conservation law. It is clear that in order to determine an exact differential on a 
particular solution surface S, a 1-form which is a generalized conservation law 
need be defined only in a sub-region of 2’ containing the image of S under the 
normal mapping. We now show that for regular variational problems, the 
generalized conservation laws can be characterized globally in terms of solutions 
of a single uniformly elliptic system. Let S be a solution surface whose normal 
directions omit the positive z direction. Let dy=@(h) be an exact differential 
on S which arises from a 1-form h on &. Let 2* be the surface w =G(g, T) 
constructed according to Lemma 6. As before, let £) be the point in which the 
positive w-axis intersects X. As noted above, in this case we need only consider 
h in the sub-region X'— f,. The projection through the origin of 2’— f, onto 2*, 
followed by the projection parallel to the w-axis onto the o, t plane, defines a 
4— 1 mapping of 2'— f, onto the o, t plane. We may therefore take a, t to be local 
coordinates in Y — fy. In terms of o, t/ is given by an expression h= A (0, t) do + 
B(o, t) dt. If V is a coordinate patch of S corresponding to a coordinate patch 
U of 5*, then in terms of the local coordinates in V, U 


(78) h = (Ao;z+ Br;)du+ (Ao;+ Br;) dv, 


(79) dy = O(h) = F*((Aoz+ Br;) dx — (Aoz+ Br) dy). 
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Let &,7 be the functions on S defined as in the proot of Theorem 3. As before, 
—,7 may be taken to be independent variables, at least locally, on S; %=@; . 
%,=t; and z(f,7) is a solution of the non-parametric variational problem ; 
6 [ G(c, t) dé dyn =0. Let p(x,y) be the function arising by integration from dy. ° 
It then follows by a calculation that y:=B, yp, =— A, and hence A;+ B, =0. 
The pair (A, B) is therefore a conservation law of the variational problem 
6 f G(o, t) dé dy =0, and hence is a solution of the system 


(80) Goya = Gon Das 2G Ave GAs a Da) 


Since the Euler equation of the above non-parametric variational problem is 
uniformly elliptic (Lemma 6), it follows that the system (80) is also uniformly 
elliptic. The above argument can be reversed to show that if (A, B) is a solution 
of the system (80), then ®(Ado+ B dr) is an exact differential on S, wherever it 
is defined. 

Let S be a solution surface whose normals make an angle =f with the positive 
z-axis, and let (A, B) be a solution of (80) in the bounded convex domain D(), . 
where D(f) is defined as in the proof of Theorem 3. The image of S on 2* under 
the normal mapping is then the region D* (6) which lies over D(A), and it follows 
that (Ado + B dr) is defined and exact on all of S. Similarly, if S is a solution 
surface whose normal directions omit only the positive z direction, and (A, 5) 
is an entire solution of (80), then ®(A do + B dt) is defined and exact on all of S. « 
Now let S be a solution surface whose normals omit no directions. Let 2* be as 
defined above. Let o’, t’, wm’ be the coordinate system obtained by the rotation — 
o'=0, t =—T, » =—w from the a, t, w system in the space Y. Let 2° be the 
surface constructed in the primed coordinate system, according to Lemma 6, in 
the same fashion in which 2* is constructed in the unprimed system. 2” can 
then be represented in the form w’=G’' (o’, t’), where G’(o’, t’) is defined for all 
o’, t’. Let £9 be the point where the positive w’-axis intersects X. As before we 
take o, t to be local coordinates in X — #,, and in a similar manner we take o’, 1’ 
to be local coordinates in X'— 4). The induced mapping a (o’, 1’), t(o’, t’) is then 
a 1—1 mapping of the punctured o’, 1’ plane onto the punctured g, t plane, which 
carries the origin in one plane to the point at infinity in the other. Let h be a 
1-form on 2’ such that @(h) is exact on S. In terms of the local coordinates a, T 
in 2'— py, h is given by an expression h=A(o, t)do +B (o, t) dt, where (A, B) 
1s an entire solution of the system (80). Similarly, in 2’— 4%, h is given by an 
expression h = A'(o’, v’) do’ + B'(a', t') dt’, where (A’, B’) is an entire solution 
of the system 
(81) Get Ag =Gye By; 2Gy er Ay = Gye (Ay + By). 


oO 


In neighborhoods of S whose normal image on ¥ contains neither po, OF fo, We 
must have 


(82) P(h) = D(A' do’ + B' dr’) = G(Ado + Badd), 


from which it follows that 


(83) A’ =(G'/G)*"(doy+Bty);  B’ =(G/G)(dop + Br,). 
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The requirement that (A’, B’) be a solution of the system (81) in a neighborhood 
of the origin thus imposes a condition on the growth of the functions A, B at 
infinity. 

In the minimal surface case, Y is the unit sphere, G(a, t) =$(o?+7?—1), 
(G'(o'", t') =3(o'?+-1'2—1), and the systems (80) and (81) reduce to the Cauchy- 
‘Riemann equations. A+7B, A’+7B’ are therefore entire analytic functions of 
@=0+1T, o’=o'+11' respectively. The mappings of 2 — p,, X — fp onto the a, t 
plane and the o’, r’ plane respectively, are obtained by stereographic projection 
through the north and south poles, respectively. The relation between 9 and 0’ 
is easily seen to be simply that 9 =1/0’. If we put /(p)=A+7B and f’ (o’)= 
A'+7B’, then the equations (83) become 


(84) fH’ (o') =|e'|*F(e(@’)) Gy = — oF (4/0’), 


where 0, is the complex conjugate of 9,,. Since f’(9’) is assumed to be regular at 
the origin, the entire function /(@) must be a polynomial of degree <2. We have 
proved the following theorem. 


Theorem 4. Let 6 { J(p) ds dt =0, be a regular parametric variational problem. 
Let &* be the surface w =G (a, t) as defined in the proof of Theorem 3. The ®-map- 
ping of the differential forms on &* into those on a solution surface S is then well 
defined. Let (A(a, t), B(o, t)) be a solution of the uniformly elliptic system (80). 
Then (a) tf the domain of (A, B) includes the bounded convex domain D(B), the 
differential form B(Ado-+Bdrt) is defined and exact on every solution surface 
whose normal vectors make an angle =f, with the positive z-axis; (b) if the domain 
of (A, B) is the entire o, t plane, B(Ado+Bdt) is defined and exact on every 
solution surface whose normal directions omit the positive z direction, and (c) tf the 
domain of (A, B) is the entire o, t plane, and (A, B) satisfies the growth condition 
(83) at infinity, ®(Ado+Bdt) can be extended to an exact differential on any 
solution surface, 1.e., in particular on solution surfaces whose normals omit no 
directions. All exact differentials on solution surfaces which arise locally from 
conservation laws of the local non-parametric variational problems can be obtained 
in this way. In the minimal surface case, the system (80) reduces to the Cauchy- 
Riemann equations, and A+7iB ts an analytic function of o+1t. The growth 
condition in part (c) implies in this case that the entire function A+71B 1s a poly- 
nomial of degree <2. 
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Some Properties of a Harmonic Function 


of Three Variables given by its Series Development 
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Introduction 
Integral operators transform analytic functions of one or several complex 
variables into solutions of linear partial differential equations. These operators 
formalize a procedure which enables us to derive various properties of the gener- 
ated functions from corresponding results for analytic functions of complex 
variables (see [6]). In the present paper we shall consider the operator 


x 
aap [AMOS X=(ny,2), 
[el=1 


which transforms an analytic function f,(w,¢) of w=x-+2 2 (Gee eeliel : 2(€—C) 


(1) H(X) = 


and € into a harmonic function of three real variables. The generated function 
H(X) is complex-valued; taking its real part, we obtain a real harmonic function 
(see (eee, by 40,09, 10, fo, 20; 17)). We assume at first that f, is of the form 


(2) fy (u, 6) = f(u) o* 


where / is a meromorphic function of one complex variable and we investigate 
relations between the coefficients of the function element of H(X) at the origin 
and the location of singularities and growth of H(X). When we attempt to 
develop the analogues of the classical theorems, various complications arise 
which have to be overcome. 

In the present paper the following questions are considered: 

I. The operator (1) generates complex functions, while we are primarily 
interested in corresponding relations for real harmonic functions, Q(X). 

In §2 we indicate how the coefficients change when we pass from the develop- 
ment of a real harmonic function Q(X) = Re[H(X)] to that of a complex harmonic 
function H(X) and vice versa. 

II. The right-hand side of (1) represents the functions H(X) only in the domain 
of association of the integral operator ((6] p. 49). By an analytic continuation H(X) 
is defined ‘‘in the large’’* (7.e., in the whole domain of its existence). When 


* As indicated in [4] the analytic continuation of the generated function can be 
ac 


pat 


‘ : : 1 
obtained by varying the integration curve ¥ in the operator Say [ Feng) 
g 
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formulating the results for the harmonic function, we have to determine the 
region where they are valid. 

III. The associate function f(w) is defined in the plane, while H(X) is a 
function of three variables. In the case when the associate (2) is rational the 
generated function H(X) is algebraic in x, y, 2; its singularities, corresponding 
to the poles of f(u), lie along circles and segments of straight lines. In some 
cases these circles degenerate to singular points. In order to formulate con- 
veniently various relations between the coefficients of the development of the 
real harmonic function, we introduce certain matrices, (1.3) * and (1.14) below. 
The coefficients of the development (1.1) of a (complex) harmonic function form 
a triangular matrix (1.13), which is a sum of one-column matrices (1.14) of the 
development 


foe) co 
H(X) = >" a, Tae (X) = ez Dd) am f umerynorde, 
n=0 n=0 [¢j=1 

6 N=mn, N=pn+x, m,p,x are fixed (real) constants. 

Using results of HADAMARD [17] and NEVANLINNA [/6], in §1 we derive 
connections between the coefficients of the development (3) and the density 
branch lines of the function H(X). See [5], p.549. In §2 we determine how 
the coefficients of the development of a real harmonic function vary if we rotate 
the coordinate system. This yields a generalization of the results of §1. In §3 
a theorem referring to the structure of certain essential singularities is proved. 
In §4 bounds for harmonic functions with circular branch lines in terms of the 
coefficients are derived. In §5 we indicate some generalizations of our results. 


§ 1. Relations between the coefficients of the 
development of a harmonic function and the density of its branch lines 
A (complex) harmonic function of three variables regular at the origin can 


be developed in a sufficiently small neighborhood ¥(O) of the origin O in a 
series of the form** (see [6] (21), p. 42) 


co 62n co an -lx—n 
r n!1' | i(x—n) « r - 

2 stra tna X) = 2 Dine Get lecnpl” Em in—nl(cosd) e' ne, XEv (0), 

n=0 x= nw—0 +0 ‘ 
where 

mwdet mt oe eee asa 
(1a) LP.) = On i | w’ 8 Sy ‘dl 
[o|=1 


are complex spherical harmonics multiplied by constants; 7, 9 and g are spherical 
coordinates. When generalizing results of the theory of analytic functions of 
one complex variable*** to the theory of harmonic functions, the following 
difficulty arises: while for functions of one variable we usually consider the 


*(1.3) =(3).0f §1,. 
** We note that i se yy 
ote that in formula (21) of [6] the summation §~ ® is used. 
axe i yA if n=0 m=—n ; 
n particular, this is the case when we consider theorems about the connections 
between a series development and the location of singularities. 


—— 
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series >) a, u”, in the case of three-dimensional harmonic functions the develop- 
n=0 


ment (1) in spherical harmonics is no longer a simple series. In this connection 
it is useful to consider at first special developments in spherical harmonics: 


(2) ; S (X) = 2g a, Ty x (X) , 

n=0 
where N=mn, K=(m+6)n+x; m,p,x are fixed integers and 0<KS2N 
for n=np, Ny sufficiently large. 


We investigate here the simplest case of series of the form (2), namely, 
ee Sig tae 
(3) S(X) = 221 ZL, nt ¢ ac. 


Let f(u) be a et i function regular at the origin O, where f has the 
development / (u =) a,u", let x be a complex number, «== a). We shall consider 


the harmonic facen 


(4) Eaies wal an 


271 f(u)—a) ° 
[o)=1 
Let A,(«), Ap(a),...,|A,(@)]S|A,41(a)| be the zeros of [/(u)—a]. By the 
theorem of MitraGc-LEFFLER, we have 


Re tial Bese sit Maw alt nae Sa 
Bel ay 2, wea FOZ oe) Hel), 


y= y= 1 


where o, are positive integers, M, are complex numbers, #,(u) are polynomials 
fo, @) 


and e(u) is an entire function. The series >) s,(u) converges uniformly in every 
7 
closed bounded set which does not contain the poles A,(«) of g,(w). Therefore, 


we can interchange the order of summation and integration in 


(6) G,(X) = Tai Lf [dst u) + e(u yor ae. 
j=? 
Whence 
oo M, (~—A %)) + [a A, («))? 4 y2+ 22)8 x 
E10.) = 2 Gres Se = [3 
ak ) Difiee (x))?-+- y2-+-22 ty+z 1 


oom 1 


2. 24 . a 

ax” | ((4—A, (a)? y?-+2?)* 
= 2 g28 

S | —(x—A,(«))+[(#—A,(a))?+y? +27] | =< P(x)} + BX) 


ai Sete ae DAs 


ty+zZ 


act 5) S,(X) + E(X), 


where P(X) is a harmonic polynomial of X=(x, y, z) and E(X) is an entire 
harmonic function. 
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We shall show that the series y S,(X) converges uniformly in every sub- 
y=1 

domain of x2+ y2+22<co which does not contain the singularities of S, (Xe: 

Let %p=(x, y, 2|2+y2+22< R%). For every R<oo, R+|A,(a)|, | 


Ute eaves © 


lo, @) 
we divide the sum )’ --- into two parts 
p= 


N lee) 
y=1 y=N+1 
so that |A,(a)|<R for y<N and |A,(«)| > for y>N. If we apply the operator 
Ba Say Eis es to the second sum, we obtain a harmonic function which is regular 
221 G 
;C|=1 


in %, since |x+¢y cost+7zsin¢/?< 27+ y?+ 22S R* for (x, y, 2)€ Xp, while 
|A,(«)|>R, for v>N. 

From (7) we deduce that S,(X) has a branch line of the second order along 
the circle 
(8) #=ReA,g), 92-2 = (Im Ages 


if ImA,(«)=+-0. If ImA,(«)=0, then the branch line shrinks to a point. We 
call this point a degenerate branch line*. 

If ImA,=0 and x=0 [see (4)], then S,(X) is single-valued and infinite at 
the point x=A,(«), y=z=O. If x0, then the function S,(X) becomes infinite 
along the segment [A,< x00, y=z=0] if A,>0, and along [—wXxSA,, 
Wesel coe A) OCR | Ol Daca ott 

Thus we see that the branch lines and the degenerate branch lines of G, (X) 
lie on the spheres 


(9) w2t y2 2 =A. (a) 
We observe that the poles A,(«) of the function g, (u)= 


= i ere 


Faye the «-points 
of f(w). Applying a theorem of NEVANLINNA ([J6], p. 72), we obtain 


Lemma 1.1. Let 
(10) G,(X)=s, [ 7 


2ni J [f(u)—a)’ 
t]=1 


where f(u) 1s a meromorphic function of u, f(0)+a, and let A,(a) be the zeros of 
[f(u) —a]. Then 


(1°) the harmonic function G,(X) has branch lines** which lie on the spheres 
(11) Cale wll Va a ee OO (e| P 
(2°) If the series 


(12) ‘ . ' A290; 
pa |.A, («)|? 


converges for three different values of «, SAY 4=4,, U=1, 2, 3, then tt must converge 
for all «, xf (0). 


- If the function S,(X) is continued to complex values of x, y, z, it has a branch 
manifold @ and A,(«)€@Z. 


** Tf ImA; (x)=0, the corresponding branch line degenerates to a point. 
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In order to be able to formulate some of the results without referring to the 
associate functions, we shall connect the study of singularities with properties 
dt coefficients a, of the development (1) of H(X). The coefficients {a,,,,} form 
1 triangular matrix 


NA 


{90 
No Ni Ye 
(13) [Gn] = 499 421 G22 Gog Gog 


439 431 432 M33 G34 A35 Ase 


which can be written as a sum of one-column matrices 


Ago 0 6) 6) 
“1 Ao Ye 0 
ey bal Were) femal Weer bua Weert (las 
; on 
ae 
(14) = [ay,u;N=n, M=n]+ [ay yi N=n—1,M=n—2)4+ 


tf @xems N= — 1, Mn) pb -- 


co 

where ay y=0 if N<0O or M<oO. In order that >’ a,u” be a function element 
. n=0 

of a meromorphic function for |~|<oo, the necessary and sufficient condition 


is that 


. l; a Ot i 
(15) line) dy em DPI 2,, 
Joo bj_-y4 n—> 0o 
ay, 7 ant» antp ay, 
2APTDy i 
Die)—|en41 Ante +++ Mtpti] = (es Ant pty sss Anti 
2 eee . 
ant p Ant p+ pres antep An+2p ve; * Antp 
p(p+1) 
SoA il) eee PRR ( (ercratrst) 1) Paar yet Rarer 
v=0,1,...,p 


where (()) is a Toeplitz matrix. (See [8], p. 335, [11].) 


Theorem 1.1. Let the constants {a,} of the development (2), N=n, K=n, 
satisfy the condition (15)*, and let 


16) F(X) = La, Ds) 
96 one-column series. Further, let 


47) G(X) FER) =L2,0) Tm X)+EX), XEKO), 


ce . . . . 
* Ya, u” is the function element at the origin of a meromorphic function. 
n=0 
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where N,=n—c,, My=n—cy. Here c, and cy are constants (integers) and 


00) = Taya) 
(18) Biju) oe [deen jaa) ce ee et 
(4) —%) RaA, 2, Ades 
az,=0 for k<0 
Gy. (ap te - O. Pors i 
(—1)¥ Ay fey Wile oon 
Coen Aay-1 4-2 4m—z --. U—% 
Ay Ay-1 4y-2--- & 
Here (( )) is a Toeplitz matrix, 
O- We Wenge we 
OP A) ae ee 
ea wee es 
O20 Ss One trearl 
Oh 2 Onsse 18, 


and F(X) is an entire harmonic function. 

Then the only branch lines (or degenerate branch lines) of the second order of . 
the function G,(X) are circles. These branch lines lie on the spheres [x?+ y?+22= 
|A,(«)|?], »=1, 2,.... The radit |A,(a)|, v=1, 2,..., have the property indicated 
under 2°, p. 210 (see Lemma 1). 

Theorem 1 can be generalized in various directions. Instead of assuming that 
we have one column, G,(X), we can consider a sum of finitely many such one- — 
column series. 


Further, one can use the associate functions 


le.2) 
(19) > a, (a Mae a 

n=0 
where N, P, x are fixed constants and |P|<N. In this case, instead of circles, 
we obtain, in general, certain algebraic curves as singular lines. 


§ 2. Harmonic functions with circular branch lines in the parallel planes 
Oye de Y ap og e= By, vied, 2535.4 
In §1 complex harmonic functions have been considered. One can easily 


formulate these results in terms of the real part of the development (1.1) and (1.2)§ 
We have 


1 no 1 ne—ntn al 
guaat ax(X) seme | Ot ey 
é|=1 
(1) tpl ay i 
(n+|x—n|)t” laa (cos 8) [cos(x — n) p +4sin(— x +m) 9g] 


def : 
“a Qn, 2[x—n| (X) +1 sgn (x va n) On a leanleieas sgn 0 = 0. 
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Let Gin» =O, 1,.0., 20, n=0, 1,2,... be real numbers, Then 
2” 
2 ine On u(X) = Re {4n0 Qno(X) a 


(2) bie a 2) ae a sgn. (x Lae n) In, 2 |x| — -1) (Oxo ra (X) +-isgn (x xy n) Qn,2a)-1(X)) 


where 
1 Aare : : 
(3) RN Nee ‘ = Gn, 2|n| A 7 sgn (— x) Gn, 2\n|—1> x= 0. 


(3) shows how the coefficients change if one passes from the series in complex 
harmonic functions I[),(X) to the development in real (modified) spherical 
harmonics Q,,,(X) and vice versa. 

The branch lines in Lemma 1.1 le in the parallel planes x=c,, where c, 
are constants. We shall derive conditions on the coefficients g,,, in order that 
the branch lines he in finitely many sets of parallel planes 

Bieta yp og t= ey, 91,2305 0=1,2,...,P) Poo. 


Every rotation of three-dimensional space, 


1) aE 
x Cains One 
y=1 


(see [74], p. 103), transforms a real homogeneous polynomial into a homogeneous 
polynomial of the same degree. Therefore, 


2n 


(5) Onu(X) =D ae, Ony(X™), 
u=0 

where d‘) are real constants (see (2)). Let 

(6) Diy = (Gin Inno eeer Gag) 

and 

(7) Gin = (Qno> “ bipees onl pee ea) 


be vectors in (27+ 1)-dimensional space. Then 


2n 2n 2n 


(8) 2) In >» dv) Q RO es Pa (qd, di, .) On; (XD), 
a=0 u=0 
where (q,,° d,,,.) —S Qn», 4%) is the inner product of the vectors q, and d,,,. Let 
(9) h(X)=2 > de Qn (X XEV(0), 
n=0 x=—0 
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be the development of a real harmonic function h(X), regular at the origin 
(see [19], p. 393). If we rotate the coordinate system (see (4)), for the transformed 


function we obtain the development 
(10) h(X) = h(X™) = D aC dys) Opp) el meee IO) 
n=0 w=0 
By (2) and (3), 
(14) A(X) =A0(x) Ref SY) a, Tn n(XO)}, 
h=0 x=—n 
where 
(4 2) He u =(q,°d nm “2 |al) St sgn fe x) (Vi 2 ||—-1) s 


Thus we see that 
wo (XK) — Ret Re{ >) 2) 


x Py tnx (COS9,) (cos (% — 2) py, + isin (x—n) n)|} : 


gin—Hl yy) 


(1) ease bse Sts 
Onn (n+|n—x]|)! - 


(13) 


Here 7,=7,9,, Y, are polar coordinates in x), y®, 2-space. We shall show that 
the series in brackets {...} converges uniformly and absolutely in a sufficiently 
small neighborhood “(O) of the origin. 4“ (X®) is a real harmonic function 
regular in (0). According to [18] and [19], p. 399, hY(X™) can be developed 
in a uniformly and absolutely convergent series 


= oy Rete Wa) a! 
AY (XM) = pes ‘| » naz) 14 Py jn —z (COS Hy) cos (x — n) gy — 


ee (n+|n—x|)! 
(14) a 
“ ae an, n—x n! n : 
= Sel Z| Y Fe eog (cos #,) sin (x — n) n). 


If we set g,=0, we see that 


15) oy 


converges uniformly and absolutely in (QO). If we differentiate (14) with respect 
to g, and then set y,=0, we see that the series 


es) de Im (ap, 4 INH) 5 | 
(16) aN ee P,, inn (C08 9))| 


Re (ays t ah 


(n+|n—x])! 


ny Py |n—x| (cos 6,)| 


converges uniformly in ¥ (0). 
In this way we obtain the complex harmonic function 


(17) HO(XOM—Y Ya Ae) Xe 70): 


n=0 x=—n 


AY (XM) = RefH®(X@)], 


where the coefficients al), _, are given by (1 De 


n,n—Hx 
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Corollary. Let q,,,, be the coefficients of the development (9) of a real harmonic 
function which is regular at the origin, and let every q,, be a sum of s numbers, 


Ss 


(o) 
(18) =e 
o=—L 


The ee possess the following property: Let 


(0,%) (2%) 
(19) ( L, i; ears +>0 for wo 
where 
(0, %o) (a, #0) = 
(20) L, = lim Dd 
n—> CO 
(0) (0) (0) 
an, ntx9 Aantint1+x%p Sieke An+pn+p+He 
; (0) () (2) 
(2, %g) Antintltyy I+2,n424xp ves Ant uti ntytitxe 

(21) Di,= e ; 

(0) (0) (0) 


Antpyn+irKy Antu+intutlt+x tes An+2u,n4+2y+Hp 


shat ( ; 
The coefficients Le and ee are connected by the relation (12). Then 


(22) h(X) = > h, (X); 


here the harmonic function 
co 8 60an 
(23) h,(X) poe X dn Q ral (X), X€V (0), Orth Dieters 55 


has pole-like singularities with branch lines (degenerate branch lines) lying in the 


planes 
(0) 


(24) Oy ¥ +n, y + Oy, 2 = const., 0 SAS 2 A .61558 

Therefore, Lemma 1.1 and Theorem 1.1 can be formulated in terms of veal 
harmonic functions. 

When considering harmonic functions of two variables, one can introduce 
a composition which corresponds to multiplication of analytic functions. More 
precisely, if H,(&, 7)=Ref,(¢), €=&+7%, k=1, 2, one defines the composition 
* by 
(25) H, (§, 7) * H2(&,) = RelA (C) f2(C)]. 


An analogous composition can be defined for harmonic functions of three 
variables. 

In (4.1a) the harmonic polynomials [,,(X), X=(%, y, 2), were introduced. 
Every real harmonic function H(X) regular at the origin can be written in the 


form 
(26) H(X)=Re{ > YAM}, XEVO), 
n=0 x=0 
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where A,,, are constants, and ¥(Q) is a sufficiently small neighborhood of the 
origin. We define the composition * of harmonic polynomials of three variables | 


by the relation 
(27) A” Lee (X) ce Ase 1 Fame (X) = Ay Ane Lite: Hy + #2 (X) *. 


Ny Hy 
The right-hand side of (26) converges uniformly and absolutely; the convergence 
of the series for H,(X)*H,(X) follows from the convergence of developments 
for H,(X), k=1, 2. Consequently, (27) yields a composition for harmonic func- 
tions regular at the origin. See [1] p. 647. 
The composition * is associative, commutative and distributive. 


Remark. The same composition can be defined by (29); see below. 
If 


a at 
(38) wa —Re | fiver Z|, #=1.2.3. 
\e|=1 
(29) A(X) = HC) + Hy) = Re] 8 filo.) felt) SI. 
\é|=1 


The above result follows immediately from the representation (6), p. 40 of [6], 
of the [7 (X). 

H,(X) will be called a qg-product (quasi-product) of H,(X) and H,(X). If . 
H,(X)=1, H,(X) will be called a g-reciprocal of H, (X). 

Lemma 1.1 can be reformulated as follows: Let G, (X) (see (10)), be a g-reciprocal 
of [F°(X)—«a]. Then the functions G,(X) have the properties listed in 1° and 2°, 
p. 210. 


§ 3. A lower bound for the radius of an essential singular circle 


To obtain further properties of singularities of a one-column harmonic function 


2 


n~*n,n+x 


(1) S(X)= 4,0 n4,(X),  XEV(O) 


n= 
in terms of the a,, we consider new combinations of the coefficients @,. Les 
the associate 


(2) S(u)==)>)a, ute™, |u| <7, 7» sufficiently small, 


of S(X) be a meromorphic function in a circle of radius R, R<oo, and s (u) 
has infinitely many poles in the neighborhood of |w|—=R. According to the 
theorem of HADAMARD, 


(3) lim es egies 


Here the J, are the quantities introduced in (1.15). The harmonic function S (X) 
has an essential singularity along a circle (or at a point) on the intersection of — 
a plane x = constant and the sphere 


(4) ty RA ane Re 
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In the following we shall give a lower bound for the radius of this circle in 
terms of the a, 

LinAy ba404 nfs A,,, v=1, 2,..., are the poles of s,(u), the function S(X)= 
iy S,(X)+£(X), see . 7), has infinitely many branch lines (or degenerate branch 
lines) Se 


(5) ey? et A 


Two possibilities arise: 
| 


} 


1. The A,, go to zero for y->ov; in this case S,(X) has an essential singular 
point. 

2. The A,, for »->oo remain larger than a positive constant; in this case 
S,(X) has an essential singular circle. 

feet 
(6) U=C,w + C,w?+ cow? +.--- 


be the development of the function which maps the unit circle onto the re- 
gion Z. Here & is the domain bounded by the segments of the straight lines 
E=—o and =— 0, u=&+7n and by two arcs of the circle |u| =R; each of these 
arcs intersects the real axis. The end-points of the arcs are 


(7) Cee (hae 02) ef Oe 

met C,, n=0, 1, 2, ..., denote a sequence of vectors 

(8) Cy = (1), C, = (41), Aes ca) CoalGy, eg peg sen), Oe 
Further let A,, »=0,1,2,..., be another sequence of vectors 


(9) Ay=(4), Ay=(@%), A,=(G,4),..., Ay = (Gq, Ge, --+,4)- 


Theorem 3.1. Suppose that the one-column development (1) of the harmonic 


and therefore S(X), has an 


function S(X) ts such that the limit 


AED 
unt Ro Pn 
essential singularity lying on the sphere (4). If the quantities 1, formed from 
a, =(A,,-C,,) possess the property that 
(10) lim (,/l,1)=1 for po, 
then the essential singularity ts a circle of radius A,=o 

Proof. By the transformation (5) the rectangle # (in the w-plane) is trans- 
formed into the circle |w| <1. The series s(w) (the associate of S(X)) is trans- 
formed into the series 
(11) s(u(w)) = > (A,- C,,) w”. 

n=0 
If the function s(w) is meromorphic in &, the function s(w(w)) must be mero- 
morphic in |w|<1. If the imaginary parts A, of the coordinates of the poles 
of s(uv) are equal to or larger than g, and a limit of a subsequence of the A,», 
=4 7624.2 equals op; then eae Gi, ae 
p00 


where ‘A are the quantities introduced in (1.15) with a, =(A,, - C,,). 
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§ 4. An upper bound for the growth of a harmonic function 
with a meromorphic associate 
The theory of meromorphic functions of one variable yields bounds for the } ' 
absolute values of analytic functions. Results of this kind can be exploited in 
the theory of harmonic functions of three variables. 
Suppose that 


_ &(U) x 
(1) s(u) = Flu) f 


is the associate of a one-column harmonic function S(X) (see (1.3) and (2)), 
where g(z) and /(«) are entire functions of order <9. 

Let A,=A,,+7A4,, be the zeros of {(w). They are ordered in such a way 
that |A,| <|A,,,|. If |A,|=|4,,,|, then Argd,< Argd 


(2) S(X) =5t, [ sw) u) Cede 
|e]=1 


is defined on a Riemann manifold-@ with branch lines* of the second order along © 
(3) = Aga, y2 + 22 = A,,, aad Dean 


As is shown in §1, s(w) can be written in the form (1.5) (with g,(w) replaced 
by s(u)), and S(X) can be obtained by interchanging the order of summation © 
and integration in 


, 1 ee 
(4) S(X) = ae if > SS) +b e(u]o ace 
[oJ=1 6 e=1 w=l 4 
According to (7) of [6], p. 47, 
Abeo(ac ace 1 —(x— A$)? + ((x—Ay)2+ y2+.22)8 
(5) = = 1 | : : = 
aes Uu—A, [(v—A,)2+ y2+ 22]? | LY+e 


We choose for A,,>0 the branch of (5) such that in the neighborhood of 
x=y=2=0 


(6) [(« — A,)2-+ y? + 22]! = (x — A,) 1 ++ ar et | . 
and for A,,<0 the branch 
(7) [(* — A,)? + y?+ 22]# = — (x — A,) [ + S aan l 3 ; 


(See footnote below.) We assume that 4,-+0, y=—1,2,.... Then S(X) is 
regular at the origin in the first sheet of the Riemann manifold. @. The domain 


of association of the representation (4) is the first sheet F, of W; F, is the whole 
space less the segments 


(8) t= Ay, Vat ee AS, ee le RP 


* If A,,=0, the branch line shrinks to a point and the function is single-valued 
in the neighborhood Of C4 OO). 
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(Those A, for which A,,=0 have to be omitted.) Let Y= A+ F be the part 
of the sheet A of W where 


(9a) A= 2 [4a + |A |e sxe Ana |Apal Ya 2 <0]; 


(9b) Ta) [A,y— [Aes Es Ass a Ayal? Ny (Ap | oe (4245) h, 
v=1 


Here # is an arbitrary number, o<p<ov, where @ is the maximum of the orders 
of the functions g(u) and f(u). 


Theorem 4.1. Let S(X) be a harmonic function with an associate ch ‘ae 


where g(u) and f(u) ave entire functions of order 0. Then in the subdomain SL 
of the domain of association of (4) and for sufficiently large r the inequality 


(10) Se Gree VC Cer 


holds. C, and € are constants, e>0, arbitrarily small. 
1 
F(w) 
1. Since g(w) is an entire function of order <9, in the rectangle |&| <r, |y| <r 
we have 


(11) log|g(€+in)|<2tt men <Ftt pew, Real. 


| 


Proof. We derive upper bounds for g(u) and separately. 


In accordance with a theorem in the theory of entire functions ([7], p. 268), 
since {(0) =-0, 

Ot é 
(12) | (u)| > |#(0)|re-” 


if |u|<oo and |w—A,|=|A,|~*, v=14, 2,..., p>o. If ¢ varies in the interval 
[O<t<2a], xy +7(vo+25)? cos(p+é#) varies along the straight line 


(13) [x = %, — (yo t+ %)? S (y? + 2*)) < (yo + 2). 
If the point (%9, Vo, 20) €:Y, then the arc (12) lies in %, and 
(14) |%o + 4(vo + 29)' cos (p + #)| S|%o + #(¥o + 20)4] 
Therefore, for r= (x?-++ y?+ 22)! sufficiently large, 
an ; ' (Fa )aere 
(15) [SCols qt f |e eset) fol aes gee 
an 5 f(w+2(y?+2?)? cos (y+?)) LOeet 
where ¢’= os +e can be made arbitrarily small for sufficiently large r. 
og r 


§ 5. Generalizations 
The consideration of §§ 1—4 can be generalized in various directions. 
I. Let H(X) be a harmonic function given by 
S 1s m4) pe 
(1) 3 dian |) Ene. mp-tm) ronan te (X) = f dale + 1)] C Lats 


— Dai 
n=0 k=0 i ; [¢]=1 n=0 
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Here p, m and x are integers, P>0, m>0. If the coefficients a, satisfy the 
conditions 1/n 


< 


1 


1j/n 1 i 
| Al? i 


Tae 


ay-1 ay 
a a 


. a 
inegoy || 


n—> CO 


Vy 1 n n+1 


(see [7], p. 323ff.), then H(X) is regular in the sphere x2+ y2+ 22<|Al? and 
has the singularity given by - 
1 siete 


2 | =e 
(2) 21 J yP(t™+1)—A 


S 
|5|=0 
Our previous considerations can be generalized to this and to analogous cases. 
In this way we obtain branch lines different from those considered previously. 
2 2 22 , : 
II. Instead of A,G=0, A,= - re "_, we consider the equation 


x2 ' Qy? oz2’ 


(3) A,®+A(r?)X VO+C(r?)G=0, GO= P(X), 
where A(r?) and C(r?) are entire functions of v?. As has been shown previously 
([6], p. 68), the operator 
1 
G(X) +f B(??, 0?) G(o?X) do = D(X), XC A(p,, ®), 
0 


Ps(G(X)) = 

(4) 7 
B (7, 62) =o B*(r?, 6?) = 0? >) B™ (r?) (1 — o?)”73, 
n=1 

transforms harmonic functions G(X) of x, y, z into solutions ®(X) of (3). The | 
domain ¥(p;,®) of validity of representation (1) is called the domain of 
association of the representation (4). The function B*(r?, o?) is an entire function _ 
of 7? for O<o0<1. It has the property that | 


(5) BY, a7) =0. | 
See [6], (10), (11), p. 68; (3), p. 64 and (11), (13), p.65. We can continue ®(X) 
and G(X) to complex values of x, y, z and introduce the variables 

(6) ae ee ee 

We write G(X, Z,Z*)=G(x, y, 2), ®(X, Z, Z*)=@(x, y, 2). From (5) it follows 
that 

(7) @(2(Z2Z*)3, Z,Z*) =G (2(ZZ*)3, Z, Z*), 


i.e., that ® and G coincide in the branch X =2 (ZZ*)* of the characteristic space. 
Instead of (modified) special harmonics used for the Laplace equation, we 
introduce the solutions 

-\ def : 
(8) ®, .(X) =p,(I,,(X)), = 0,4,2,..., %=0,1,..., 20. 
Ina sufficiently small sphere (with the center at the origin) which lies in the 
domain of association /(p;, ®), @(X) can be represented in the form (4), since 
G(X) can be developed in series of spherical harmonics. Therefore, 


iS) 


0) HX) = Say®(%), KEV). 


The series converges uniformly and absolutely in the sufficiently small neigh- 
borhood W(O) of the origin. It follows from the representation (4) that in 
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A%(p;,) the function @(X) is regular in every simply connected domain which 
‘contains the origin and in which G(X) is regular. The locations of the singu- 
larities of ®(X) and G(X) are closely connected. In particular, if G(X) has a 
branch line, in many instances ®(X) must be singular along the same line. This 
shows that many of the theorems about the connection between the coefficients 
of the development (1.1) and the location of the singularities of the harmonic 
function G(X) yield theorems about the connections between the coefficients a, , 
of (9) and singularities of @(X). 

As an example of a singularity of (3) which one can obtain using the integral 
operator p,;, we consider the solution ®(X) 


ps ([(~ — At)? + y? + 22]-4) 


=[(x —Ai)?+ y24 ay [ B*(r?, 0%) 62a Patches eh 
0 


(10) 


[(wo2— A 1)24 y24 2218 
If we approach 
(11) j= |e 0, yi Ae 
along *=0, y?+2?< A?, the function [(x—A7z)?+ y?+22]-4 goes to oo. 
We shall show that the second term of the right-hand side of (10) remains 


ee oeny, bounded if we approach 4 along the same path. 
1 


(12) ! AZ ey odo 1 i) ; B*(r2, t) dt 
[ (40? ee J 42 G4 22 ot}? cs ¥[(t™—T(¥, %)) (T-T, (*, n))} ; 


where 


wlun)=six—(— +H, non)=4 feat (-— eth, 
(13) 
(<4 rar — 5 +, 7 = (y2+2)8>0. 


For x=0, we have 


1 IBS A Gea) dt 1 eda * 
iy 2 | A A\t| = 2 A |t* 
Vy (7 a (=+ >) F Pie = 


1 1 


Hor 0O=7=1 and 0<A == Hee A, we have 
4 2 


r+ = 7 >> 0; Make Vac! t)| SIBD< 68, 
1 


Als 
im ar 
1 
Poinded when we approach 4 along x=0, y?+22< A?. 


In the corollary on p. 215 the assumption is made that the coefficients q,,, 


Since | ee < oo, it follows from (14) that the second term of (10) remains 


(2) (Q) : 
of the development (2.9) are sums of s constants 4, and that the das satisfy 
as limit relations. (See (2.18)—(2.21).) Suppose the series development 


ee . . 

> s Inu Dysst ) is given. Then arises the problem of determining whether the 
n=0 u=0 

Yn, can be represented in the form (2.18), where the iy satisfy some conditions 


of HADAMARD’s type. The author hopes to return to this question in the future. 


This paper was prepared under contract Nonr-225(11) and NSF Grant 10375. 
The author thanks Dr. SictaK for his help in preparing the paper. 
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Al Note on Harmonic Functions in (p+ 2) Variables 


R. P. GILBERT 


Communicated by A. ERDELYI 


| 
| 


The purpose of this note is to present two integral operators which transform 
functions of (p+) complex variables into harmonic functions in (p+ 2) variables. 
Both of these operators are generalizations of the WHITTAKER-BERGMAN operator 
B;(f, &, X°), which transforms functions of two complex variables into harmonic 
| functions in three variables [J], [2]. 
The first operator, an immediate extension of the Bergman operator, has 
the form 


o= 44 Ze (6,4 a (6 =) (6s a ical 


+ SE seal a) 6g) ba galt g) 


PG t, 
$+ Lo pa(s— B)(G— ne (Cp rae malt 
p 


the domain of integration “= iii is a product of regular contours 4%, in 


the ¢,-planes, |X —X°|<e, X= oe Lots i teen =a Spree Heals coma 
€> 0 is sufficiently small. 


o has been chosen in such a manner that o” is automatically a hyperspherical 
p+2 

° o2 

harmonic, and consequently one has » : 

rest Ox? 

the operator Bs.9(F; £; X°) transforms the functions F(o;¢) into harmonic 

functions H(X) if we consider o” as the generating function for a class of harmonic 

polynomials. We introduce these homogeneous, harmonic polynomials of degree 

+n 
n, h,(X;m), by the relation** o”= >) h,(X; m)¢”. With the residue theorem 


m=—n 


F(o;C)=0. We may realize how 


* Here and in future instances we shall use the contracted notations F(c; ¢) 


- ac. - ; 
and a to represent (0; 7,,C.,...,)) and ae 22). Ee , respectively. 
G Cy Ge dlp 
+n 
*x Here we have introduced the notation ») h,(X;m)¢”, which is to mean the 
+n +n +n m=—n 
p-tuple sum > Dy, S00 DE fy (X; m,) Cy" Cy... Ch”. This abreviation will be used 


m=—N Ms nN Mp 
throughout the remainder of Re work. 
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we may then show that the functions of (p+1) complex variables 


map onto the harmonic functions 


(oe) +n 


H(X) =Biyo(F; £;0)= D> a,(—m)h, (X; m). 
n=0 m=—nN 
There is, however, a major disadvantage in using this operator, in that o” does 
not generate linearly independent polynomials*. Because of this we now turn 
to another integral operator, which does generate linearly independent harmonic 
polynomials. 


We introduce the operator 


BF; £:X)=4 f FO) |2( ales (gis cok Ok 153, 


Dyes C 
L 
= 1 see ) 7 
T —| sm RENO 
k 4 he Ch ie O,, ( NS ? iP bie 
- Se cee ae be 
T=hr, and t= ae (C, z,) Sin Bp 
where 
i 1 
1, = cos, + sin (C+), ih A ese, 
|X — X°||< Bs 9 GAC rie Pare 4) XS (24, My s+.) Xora), 


é> 0 is sufficiently small, and the domain of integration 4 = I Y,, is a product 


of regular contours Y, joining +1 to —1 in the €,- ninkees The Or, Y, and 7 
are hyperspherical polar coordinates, and may be defined by [3] 


Lies tiCOS ie 
Ko =7 sin B, cos, 


X%3 =rsin #, sin J, cos Bs, 


Xpi2=rsin sind, ...sind, sin ep, 


where OS OS (PAO ap amet yS2a,; the constant A is defined by 
p-1 p-1 
A= Ie A,4= = at ate = 
8 k=0 (mz—mp+4)! P(mintS—2)| 


* There are only (2n+) Aer ea 
of degree m [3]. pin! 


linearly independent harmonic polynomials 


Harmonic Functions of Second Variables 225 


This operator was developed by considering the following functions, which 
form a complete, linearly independent set of harmonic polynomials of degree 1 [3]: 


| ; : __ [{ *pt41 » Xpto\—Mp 
: H(m,; +; X) ( ry +1 au Dre 1 Orr (%x/%%—1) 


1 k 
(eux = 4 - | 5 ae and Y= My_1— ™,), 


iy 
= Y(m,; O,; p) = 7" ct? [FE (sin 3,)”* Cot (cos ¥), 
: k=1 
where 


LSE 2 2 yh re 
1, = (Magi Meg2 +o + Xp y2)', 1 = 19, 


k=0,1,...,p,andn=m=m,=---=m,=0, and where the C; (cos #,) are Gegen- 
bauer polynomials. A harmonic function in (p+ 2) variables which is regular 
in a neighborhood of the origin may be expanded in terms of these polynomials 
as a series * 


m 


aS 2B, (m) 7" Y(m; 9; + 9) 


n=0 m=0 


that converges within the largest hypersphere inside of which there is no singu- 
larity of H(X). The operator B}..(F; ¥, X°) is then obtained by expressing 
the H(m,; -+-; X) as combinations of the integral representations ** 


C4 (cos 8) = 


ae AL ete Pf eos0 isin d cost iGing)t4-2d¢. 


An inverse operator for BS .5(F; &; X°) may be obtained by making use of 
the orthogonality property of the H(m,; +; X) on the unit hypersphere [3]; 
that is, 


fe 2 
Sf A(m,; +; X/r) H(m; +; X/r) = (22) iby One Ey, (Mp—1>™) » 
Q 


k=l 
where 


Ore PT (mp rtMmatp+1—k) ; 
| ee Aa 
[rm ep E ao (70, — mz)! ir r(m.+ P = ) 


Ey (My,-1, ™%) = 


ie Bass 
Pgh ee ey 


For instance, let H(X)=B5..(F; £; X°), where F(z) is of the form 


2 


m 
EG@)= >, D)-0, (1) (n3.2); 
n=0 m=0 


mv’ 
x The expression >) B,,(m) 7" Y(m; 6; +) is meant to designate the p-tuple sum 
My Mp—y m=0 
3 D+ By (1g) ¥ (rg; 8; = ¢).- 
=0 m,=0 mp=0 


x* It may be seen that the H(m,; +; X) are generated by Boi2(F; £; X°) from 


P 
the terms p(m;t)= J] t* 
k=0 
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then there exists an inverse operator 


1 Ss en gS 
BB a (F; 2; X94 = Cy (Hi; Y) =F) = [[ HO KCTS £9) dQ. 
J | 


(where the integration is over the unit hypersphere), which maps the harmonic 
functions in (p+ 2) variables H(X) onto the analytic functions of (p+ 1) variables 
F(t). This may be seen by choosing 


Tole oe hs y X et”? (t sin 3)” w CY (cos 8) | 
n=0 m=0 
(=m z ; s Y,=™,_1—™,, and k=1,2,...,9), 


applying the orthogonality properties, and noticing that 


nA 
E, (mz _1, Mx) = - p hes ie 
eae cay 


We shall now state without proof a theorem concerning the singularities of ~ 
harmonic functions generated by an integral operator Bi. .(F; ¥; X°), (j=1, 2). 
The proof is essentially the same in structure to those developed by the author 
for three and four dimensional harmonic functions [4], [4]. 


Theorem 1. If Z?*?=E{S(X;f)=0} is the singularity manifold of the integral 
operator Bi, (F; £; X°) (j=4, 2), then the harmonic function H(X) = Bt a(F; L; Xe 
is regular at the point X providing this point does not lie on the intersection 


E{S(X;0) =0}0 El 4 ¥ tae =o}, 


ot, c= OfR at, 
where the (,=wy, (C,) ave arbitrary functions of C,. 
We note, that certain special cases occur when the C, (k=2, 3,..., p) are 


independent of ¢,, for instance a particular class of singularities may occur 
when X lies simultaneously on the surfaces 
Clee ty 5 S(X0) Shc See 
v 
We now turn to a few examples to illustrate the use of Theorem1. For 
instance, let us consider the case where F(t) is a function of t, alone, that is 
F(t) =f(t»), and where f(t) has a singularity at ty =a. By making use of the 


identities COS 0, = Xgiy,_,, SN O,=7p),,_, (R=1, 2,...,p), We May represent the 
singularity manifold as 


S(X34) = (%— Do =t mt 28 41) =0. 


Eliminating ¢, between S(X;¢,)=0 and = 0 yields (x,—«)?+73=0, which 
1 


was to be expected, since these are the singularities of the generalized axially — 
symmetric potentials [6]. 


As another example, let us consider the case where F (t) is merely a function 


ba 1 -2 . — 
of 1,6) , (1S7Sp—1), and has a singularity at (G—] ee Here, 
7 
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he singularity manifold of the integrand may be represented by 


S(X; Car Gi ea) = Cj (Cia Moat s V541 [Cf Stik 1)) =i 
+ Aaljaalby 29+ $7 [E7 + 1]) =o. 


he only possible singularities of H(X) occur then for X simultaneously on 


(X;6;, «(E;)) =0, and aS (X;¢;, «(C;)) =0, where ¢,.;=2(C,) is an arbitrary 
o] 


unction of ¢;. For the purpose of illustration let = (¢;) =¢,; then upon eliminating 


arom S=0, ion one obtains the locus of possible singularities 
7 
(% 41 + 4a x) + (V4.4 a 4a7;)? = 5 

In an earlier paper [6], the author showed that there are certain distinguished 
ints associated with an operator which transforms analytic functions into 
eneralized axially symmetric potentials; indeed, it was shown that a branch 
f a potential function might have singularities on its axis of symmetry, which 
lo not correspond to singularities of the analytic function*. This situation may 


e seen to occur here also for the case where S(X;¢;,¢;.,) is given as above, 
os 


nd we consider the class of possible singularities for which X lies on S=0, ae 0, 
; 7 
nd = =0. One then obtains 4) = %j42= ++: = %i2=0 (7 =1,2,...,p). 


Ve note that this is a close analogy to the axially symmetric case, and it occurs 
ndoubtedly for the same reason. 
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Application of Conformal Mapping to Boundary 
Perturbation Problems for the Membrane Equation 


LEE A. SEGEL 


Communicated by C. C. LIN 


1. Introduction 


Lin has suggested [5] that conformal transformations are useful for boundary 
perturbation problems involving general differential equations. It is required — 
that the boundary conditions be given on some curve B, near to a curve Bg 
for which the problem can be solved; the conformal mapping of By onto B, 
will then differ only slightly from the identity mapping. Since we consider 
equations more general than Laplace’s, the transformed equation will differ from 
the original one, but by a small amount, so that the solution can then be deter- _ 
mined by standard perturbation methods. 


In another paper [10], the author has shown that Lin’s suggestion leads to 
an efficient method of solving certain problems of current interest in fluid 
mechanics, for example the two-dimensional time-dependent flow of a viscous 
incompressible fluid forced by a circular cylinder vibrating about a point slightly 
displaced from the center of a surrounding circular cylinder. In this paper we 
report the results of work on the eigenvalue problem 


RAW =0cin-Ry: WO oni Bs, (1) 


where R, is a region in the z-plane (:=x-+7y) bounded by B,. This problem 
is the simplest one which includes all the features we wish to consider but is 
already sufficiently general to describe significant physical situations. For 
example, the discussion in [/] on the validity of the Born-Oppenheimer approxi- 
mation requires treating (1) when B, is a parallelogram and B, is a rectangle. 
The work which follows, particularly §5, would be at once applicable. 

Several methods have been proposed for finding eigenvalues in boundary 
perturbation problems. (See [7], Chapter 9, and the references given there.) 
Although each of these is useful in appropriate instances, to the author’s knowledge 
no convergence proof has been given, and at least one method generally soon 
diverges. In contrast, we give conditions guaranteeing convergence of the 
mapping-perturbation (MP) method for sufficiently small perturbations. For 
a large class of interesting cases, these conditions are expressed simply and 
geometrically. To take a specific example, the method definitely converges for 
the principle eigenvalue of an ellipse of eccentricity less than about 35° 
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Often more important in practice than the convergence of a perturbation 
method is the degree of approximation obtained from the first few terms. For 
this reason we compare the MP method with three other methods for attacking 
(4) when B, is an ellipse of small eccentricity. 

Finally, we show that, when B, is nearly an isosceles right triangle, the 
2igenfunction perturbation series in the original variables must have certain 
singularities which are absent for the series in mapped variables. It appears 
that perturbations of the independent variables by conformal mapping can offer 
analogous advantages for elliptic equations to those known to be possessed [5] 
by perturbations in characteristic independent variables for hyperbolic equations. 

We begin by outlining the MP method for (1) and applying it to ellipses of 
small eccentricity. 


2. Some formal calculations 
Assume that an eigenvalue 4 and the associated eigenfunction W) are 
known for the region R, “near’’ the region R, of (1). The formal procedure 
for obtaining / is then to transform the region Ry (taken in the € plane, where 
>=&+1n) by means of the conformal mapping z=z(C; «) into the region R,. 
With this, (1) becomes 
Se toe tiie=o Ing, w= 0) ON By, (2) 


where w(&, 7) =W[x(&, 7), y(&, 7)]. The Jacobian J is assumed to have a series 
expansion of the form 


| a 

= le 
so that we have transformed the original boundary perturbation problem (1) 
into a standard perturbation problem. Formulas for obtaining the series 


(5,1) + € Jo(E,m) +++, (3) 


A=/1O + oe AM +--+; (4) 
w=wO+ew)+..., (5) 
can now be derived by a straightforward extension of the usual perturbation 


calculations used on ordinary differential equations, as in [13], p. 219—223. 
As an example, let us consider the region bounded by 


zz — (=) (2+ 2) =1, (6) 


z 


% \8 date 
which is an ellipse of eccentricity (- =} . The required series (3) for ze can 


be found from formulas in [4] expressing the mapping of €?+7%<1 onto the 
region bounded by the analytic near-circle 


2€ 


z@+eF(z,z)=1, fFanalyticnear 2z7=4-. (7) 


5 . 
The mapping is oie by z=C+ : 7 eo Ais : +. =)+ --» so that, with 


-=oexp(t 9), fa "14360? Sgecbiiiee 
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Let A,,,- and 4,,,, denote the eigenvalues associated with the eigenfunctions: 
Tin (inno) cosmap and J,,(%,,e) sinme, J, being the m‘>-order Bessel function; 
with fio.) 0; 7=41, 2, %.. In therseries 


Anne - Oe Te Dee ahaa hans = eae tae ope ae ees 
application of the appropriate formulas to our example shows that 


i, = A, =O, med, 
2 2 * 2 1 4(0 (8) } 
Ains= Mae Cade (Can) Jo? A (%1,0)40=s4in, 


the integral being evaluated with the aid of [75], p. 138. 

It is natural to compare this example of the MP method with the classical] 
use of Mathieu functions, which involves the conformal mapping of the interior 
of the ellipse onto the interior of a rectangle. As in [6], p. 212, the k‘® approxi-- 
mation to the lower eigenvalues can be found by computing a root of a poly-- 
nomial of degree 2(k —1). Another natural attack is to employ the Joukowskii 
transformation, z=¢€-+-c?¢—1, with which it is possible to work out a perturbation} 
procedure for ellipses of small eccentricity. One computes the first approximations 3 
to the eigenvalues by determining the zeros of an expression involving Bessel | 
functions. Since the mapping is of the interior of an ellipse onto a circular - 
annulus, however, these equations are somewhat complicated and are useful onle ' 
for the problem of a region bounded by two confocal ellipses. 

In order to compare the Mathieu, Joukowski, and MP methods, we have | 
calculated A,,;, and /,,, for the ellipse (6) when ¢=0.0204; the results are given 
in the Table. It is seen that the MP method is an efficient one even for this : 


Table. 4,,, and A, 


Approximation 


Method 


JOUKOWSKI 15.05 _ 
MP 14.67 | 14.53, 14.83 — 
MATHIEU 14.39 | 14.53, 14.87 14.53, 14.84 


special problem: it seems to converge better than the Mathieu method, although 
the calculations at each stage for the latter are easier. In addition, convergence 
of the MP series is proved below, while the convergence of the Mathieu method 
seems never to have been proved. 

It is natural to ask how the calculations would proceed under a non-conformal 
transformation. One such approach is to write (6) as r=(1—ecos 29)? and 
then to make the change of variable 7(4— ¢ cos 28)!=0, ¥=q. The differential 
equation becomes 


in2 2 
1 — ecos2 ein eee 
| a ee? ci 
1 2€ COS 2 + e7(sin? 2m — 2)]} dw 
—|41—ecos2 P 
| ae , Fee eee (9) 


1 Cw rw 
+ ta [1 — € cos 29) a 4 2esin 29 ee LAw=0. 
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The boundary condition is w=0 on g=1. Proceeding with the calculations, we 
again obtain the MP result (8), although the expression leading to it is different. 
This is not surprising since (9) contains ¢ analytically and should therefore lead 
to convergent series for the eigenvalues. The complicated nature of (9) would 
make a proof difficult, just as it makes calculations difficult. 


3. A convergence theorem 
The simplest convergence result is the following: We assume (at first) that 
the unperturbed eigenvalue is simple, and that B, is composed of a finite number 
of curves possessing a continuously turning tangent such that the angle at which 
two such curves meet (as seen from inside R,) is non-zero. We further assume 
that the mapping function z=2z(¢, e) is such that for some constants c and K, 
the expansion 


ge| = tt enn) + teen) + (10) 
converges in the mean to oe , the coefficients satisfying 
sl GeAe os EOC Ti — le ere (11) 
If d is the distance from 65 to the reciprocal of the nearest unperturbed 
eigenvalue, and J"is defined as in (18), we define e, by 
a= (5) @+4PKy. (12) 


If | e| is less than ¢,, then (4) converges, (5) converges uniformly for (€, 7) in Ro, 
and the terms in (4) and (5) are correctly given by the formal calculations. 

As was discovered after the present work was completed, TircHMARSH recently 
discussed eigenvalue problems for perturbed differential equations (not for per- 
turbed boundaries). He proved for (2) that if the function represented by / 
satisfies certain inequalities, then the formal perturbation calculations are valid 
[13]. TITCHMARSH’S results would require modification to suit the purposes for 
which we shall need them, but his work is sufficiently general to render inadvisable 
a full presentation of the proof of our theorem. Consequently, we shall limit 
ourselves to an outline of our argument and refer any interested reader to [17] 
for details. Our proof is an application of known theorems in the perturbation 
theory for linear symmetric operators in a Hilbert space. This theory was 
originally expounded in a series of papers by RELLICH [9]; we use the refined 
version of Nacy [72]. (TrrcHMarsH’s work makes no explicit use of Hilbert 
space theory but “‘is adapted ... from the analysis of RELLICH’’.) 

In order to apply the RELLIcH-Nacy theory, instead of considering (1), we 
consider the equivalent integral equation 

W(x, 9) =ASS G.(%, 95 0,9) W(x, y') dx’ dy’, (13) 
Re 
G, is the Green’s function of V? vanishing on B,, the primes denote variables 
of integration, and R; is the (x’, y’) region corresponding to the region R, for 
16% 
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“= 


(x, y) ([3], p. 365). Changing variables of integration, and providing symmetry 
by defining a new kernel 


K(E,n3 &,n'; €-) = GolEns €') En: €) JEG ei, (14) | 
we can write (13) as 
Sf KE ns 0; ©) w(E,9'; €) dé’ dy’ =Aul§,n; @), (15) | 
Ro 


where u=w Jt and A=A7+. We define the operators A(e) and A; in A(e) =), A;é ' 
by identifying (15) with 
A(e)u=Au (16) 
and base a Hilbert space on the inner product and norm 
(u, v) ie n)v(é,n) dédn, |u|? =(u, w). 
The bound on the 4; eit by the RELLICH-NAGy theory, 


|4;|<Ky,ct, K, and ¢, constants,  |A,|| = lim sup | ( (Au, u)| (17) | 


||u|]= 


can then be deduced from (15) by estimating the integrals involved, such as 


r= ff f[/Gen; &, 7') dé dn’ dédn. (18) | 

Ry Ro t 
The remainder of the proof consists largely of showing that convergence in the : 
mean for the eigenfunction series implies their uniform convergence and that the : 
formal perturbation calculations from the differential equations yield the solutions i 
whose existence has been demonstrated. 


We now mention three generalizations of the principle result: | 

a) Even if the unperturbed eigenvalue 4 has multiplicity uw, 1<ju< ©, the: 
appropriate formal calculations can be shown to be valid for sufficiently small | e|._ 
Explicit estimates require further information on how the several perturbed 
eigenvalues diverge from 4 as |e| increases. An easily handled case is that 
considered by Nacy [12] where an unperturbed eigenvalue A® of multiplicity w 
splits at order ¢ into yw simple eigenvalues. 

b) Using an appropriate Green’s function, one can show without difficulty: 
that our results hold if, in (4), the Dirichlet boundary condition W=0 on B, is: 
replaced by either the Neumann condition de =0 on B, or, if B, is taken to be: 
composed of two arcs B” and B®), by the mixed boundary condition oom 
B®), W=0 on B®), oF 


c) One obtains a better convergence estimate than (12) by adapting the work. 
in [12] to show that if, in addition to (17), the inequality 
|A(e) — Ag|< M(e) (19) 
holds, then the series for A and w of (16) converge if | 


lel<at, M()< 4. (20) 


| 


. 
| 
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4. Some regions for which the method converges 

We shall first demonstrate that the MP method converges for boundaries B, 
which are analytic near-circles. For this class of curves we prove in the Appendix 
that an inequality of the form (11) is satisfied; this shows that the formal calcula- 
tions are appropriate for | e| sufficiently small. For simple eigenvalues, an estimate 
can be obtained from (20) with the aid of the following result of WARSCHAWSKI [14]: 
Suppose that the boundary B of a region R is a closed Jordan curve represented 
by the equation y =7(#) where 7 is a continuous function of period 2a. Suppose 
further that, with F(9) =" ae for some 6 


1<7(9) <1+6, |F@)|<6, |F(@+7)—F(d)|<|ctlo. (21) 


Then if z =2(€) maps |¢|<1 conformally onto R such that z(0) =0 and 2’(0)>0, 
the inequality 


lz" (C) — 1] S2(1 +.0)2 4? exp (6%) —2— d= (0) (22) 
holds for |¢|<1. 


Using this, we see that in the expression 
ILA (2) — Ag] uP 
/ , & d / / , , , / 2 
=SS{ff eens.) || se 80] | ge 0)| — 1] eG 0) a8 an aban, 
1D ie 


the quantity in the square brackets is less in absolute value than Mj -+-2M,=M,. 
In the Appendix it is shown that the remaining part of the integral is less than 1, 
which yields 

|\[A(e) —A,] wfP?S MP ule, |A(e)—4,|S™%, (23) 


so that the M of (19) may be identified with M,. 
As an illustration, we apply this work to the principal eigenvalue of the 


ellipse (6), which we write in the form 


E:r=(1— «cos 28). (24) 


Consider c=0. As (1+¢2)-$SrSX(1—.e)~, to use WARSCHAWSKI’s theorem we 
first subject (24) to the mapping z, = (1+ aua)hag which takes E into an ellipse EF, 
satisfying (21) if 6 =2e(1—«). Let z,=2(¢) map |¢|<1 onto the interior of 
E,; (21) then holds for 2 (¢). This leads to 


|2(0) —1|< (M41) (1+e4#-1=M. (25) 


The lowest unperturbed eigenvalue, 7.e. the lowest eigenvalue for a circular 
membrane, is 5.78, and the next lowest is 14.68, so that, in (20), d=0.105. Using 


My of (25), M,=5.86 +0(6%), 6 aa and M=(M,)?+2M,, we see that the 


requirement M< z is satisfied when e<0.0012. From the Appendix, the other 


requirement of (20) certainly holds for this range of ¢. Therefore the principal 
eigenvalue of the region bounded by the ellipse (6) is an analytic function of « 
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provided 0<¢<0.0012. Equivalently, since the eccentricity E is related to e: 
by E2 =2¢(1-+e)+, the lowest eigenvalue of the elliptic membrane is an analytic : 
function of the square of its eccentricity provided that its eccentricity is less } 


than about 545. 


We remark that, by use of results in [4], it can be shown that the MP procedure : 
is convergent for any membrane bounded by a curve departing slightly and | 
analytically (in the sense made precise in [4|) from some non-singular curve. We 
note that in any rigorous discussion it is an advantage of the MP method to be | 


able to draw on the large body of research in function theory. 


5. Triangular membranes 


Exact solutions to the membrane problem are known if the boundary is a 
rectangle or an isosceles right triangle ([7], p. 754—756), or a regular hexagon or 
an equilateral triangle [2]. It is natural to hope that these solutions, too, can be 
perturbed successfully, although one at first would perhaps be pessimistic because 
of the mappings’ branch points at corners of the polygons. We consider here the 


particular case of a triangle which is almost an isosceles right triangle, and we | 


prove that the perturbation method yields convergent series for the eigenvalues 
and eigenfunctions. It would seem that the proof could be carried through for 
any region bounded by straight lines. 


Consider a nearly isosceles right triangle 7. whose angles at z=1 and z=0 


are e +aen and ae + Pea (« and f are constants). The MP procedure requires © 


mapping 7) onto 7, where 7) is a triangle whose angles at €=0 and(=1 are 
both ae By the relation 


h if 
=O, [Ey Pere PtOede Cot = f (4 — 8) ee age ee 
0 0 


the interior of 7, is mapped onto the upper half /-plane, the points z =0, z=1 


remaining fixed. With the same fixed points, 7) is mapped onto the upper — 


h-plane by 
GaGa) EP elma de, (27) 


which completes the originally required mapping in terms of the intermediate 
complex variable h. As the eigenvalues and eigenfunctions for 7, are known, the 
MP method may be used. Integrations will be performed in the /-plane, so we 
note that if h =a+7b and F(h) =F {¢), then 


es “=| dC |2 os an Bs 
[[ Fagan = [[F “~\dadb = {fF \H ]1—h|-3dadb. —(28) 
To b=0 b=0 
The convergence theorem requires that te satisfy (11) where Z| = .Cy aaa 
ay G : 


A=|h|’|1—h|*. To demonstrate (11), we first consider the formal expansion 


At = pS gl H,, H, = (In zl ; 


n! 
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A bound on 


: a | 4a db can be obtained by dividing b=0 into region I 


b=0 
where A<1 and region II where 4=1, and by observing that for any y>0, 


4 

(In aS ea pie Oe ts (29) 
/ 4 

= aa ” og? {prs ea (30) 


Stirling’s formula, as in [15], p. 251—253, gives (m!)*=4n2(m +1)4"t2e—4"+9), 
which shows, using (30), that for some constant K 


Jf Blades Gy ffl 


Application of (29) shows that an inequality like (31) holds over region I. Choos- 
ing y small enough to insure convergence of the bounding integrals, from (28) we 
have i i HedédnsK ee Similarly, from (26), the coefficient of «” in se is 


€ 


—#+78) 4 — pl 8 dadb. (31) 


less hen a constant multiple of ek So is the coefficient of ¢” in — for 
0 


o< 7 <1 implies that the two series have the same radius of convergence. 
The desired inequality, i f gp, dédn<K (= y ", now follows without difficulty. 


We have therefore een convergence of the method for the triangle T,. 
In particular, if | ¢| is sufficiently small, any eigenfunction may be represented 
by a power series (5) which converges uniformly for all (&,7) in 7). With the 
series €(x%, y; ¢)=x+e&(x, y)+--- and n(x, y; e€) =y+emn(x, vy) +---, deriv- 
able from the mapping and the definition W(x, y; ©) =w[€(x, y; €), n(x, y; €)], 
one can find expressions for the terms in the expansion 


W(x, y; 2) =Wolx, y) + eM (%, y) +o (32) 


of the eigenfunction in the original independent variables; e.g., 


W, = wp, My, = 85) | ma (5) + Wy. (33) 


The expressions for succeeding terms in (32) contain progressively higher order 
derivatives of terms W; of progressively higher subscript, so unless all the terms 
of (5) are infinitely differentiable, (32) cannot converge uniformly. Even if it 
does converge uniformly, (31) compares unfavorably with (32): as is easily seen, 


a (say) is singular at the triangle’s corners so — from (33) — w, would not even 
a once differentiable while W, would be infinitely differentiable. 
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Appendix: Inequalities involving Fa for mappings of analytic near-circles: 


In 1946, GoLusIN proved that the power series in e giving the mapping: 
z=2(C; e) of |f|<1 onto the interior of the analytic near-circle (7) converges iff 
|e|<e*, uniformly for |¢]|<1. A lower bound for e* is obtained in GOLUSIN’s: 
proof ([4], p. 391 ff.). This bound is worked out in detail for the ellipse (6). There : 
is a slight error; the answer should be e*=0.0145. To extend GoLusin’s result t 


& , we note that, for a fixed e, z(€; «) represents a conformal] 


mapping and so is an analytic function of € if ||<4. Since a function separately 
analytic in each of several variables is analytic in all together, if 6 is any small! 
positive constant, then z can be represented by an expansion 


to series (10) for 


BG; é) Gee Bee (34) ) 


which converges uniformly for |¢|<1—, |e|<e*— 6. But since z =2z(C; e*—6)) 
maps a circle onto an analytic curve, z(¢; e*— 6) is an analytic function of ¢ for: 
|¢|<4, not just for |f|<1 ([8], p. 186). Therefore series (34) for z(¢, e*— 0) con-- 
verges for |¢|<1, (34) converges uniformly for |¢|<1 and | e|<e*—6, and for: 
dz dz | 


these values of € and ¢, we may obtain, from (34), series for ade? de’ and | 
dz |2 & Sey 3 
i= or . Assuming no further restriction is needed to insure that |J—1|<1,, 


(10) converges for |¢|<e«*— 6, uniformly in €%+7?<1. Consequently lim sup: 
1 


(§+-4?<1, noc) |p,|" =(e*— 6), so that there exists a constant } such that: 


IJ @nd& dy <d(e*) 4", 


SPs 
which is the inequality required in the convergence proof. 


Bound on a Green’s function integral. If G is the Green’s function of’ 
V? vanishing on &%+ 72=1, we give the upper bound for the quantity J =- 
IS a G°(E,n; &', n') dé’ dn’ dE dn required under (22). As in [3], p.377, if: 


gh = s8+n® and (9')?=(6')?-+(y')%, then G=—}2 (Ing,+Ing, 9’). Here @. 
and Q2 are the distances of (&, 7) from (€’, »’) and from its reflection in the unit: 
circle. It is easily seen from a figure and use of the law of cosines that- 
1— ee’ S@0'S1+00'S<2 and |e—o'|<e,<e +o’ <2. Using (In x)?< 2.3 x-} for’ 
te “S14 and Schwarz’s inequality, we find J2< (422) (2?) (2.3) [,+J.+ (1, 1.)%] 
where : 


\ eee 
=2f [°° aed M3 ().702 and ={f eo! re 
J (e—e') : : Jd — 00’) dg do’ 0.374, 


so that J2<0.99 and I< 1. | 
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On Lancaster's Decomposition 


of a Matrix Differential Operator 


A. M. OSTROWSKI 


4. In the theory of linear vibrations of damped systems a matrix differential | 


equation of the form 
1 (4 4, BL AClE HOO 
” dt? dt se ; 


with (1 <n) matrices A, B, C and (n x1) vectors &(#),¢(¢), usually plays a central . 
role. Closely connected with this equation are the non-linear eigenvalue problems 


(2) (AR+BALC)E=0, n(A#+BA+C)=0, 


where é,7 are the eigenvectors and / is the corresponding eigenvalue. 

In this connection P. LANCASTER [1] has recently discovered a decomposition 
of the operator in (1) which holds under fairly general conditions and appears 
to be very useful. 

Since LANCASTER’S proof is imbedded in a rather cumbersome general invest1- 
gation and is not very accessible, I shall give an independent formulation in 
Section 2 and a very simple proof of this result in Section 3. 

The decomposition in question uses only one half of the right and left eigen- 
vectors of (2). In Sections 4,5 of this note I give simple expressions for the 
remaining eigenvectors and use these expressions in Section 6 to present LAN- 
CASTER’S solution of (1) by quadratures in an alternate form which may in 
some cases be easier to handle than the expression given by LANCASTER. 

It has been found convenient to use notations different from LANCASTER’S. 


2. Theorem. Consider the differential operator 
(3) AD? AB DrsC: 


where A, B,C are constant (nXn)-matrices and D is the differentiation operator 
with respect to t. Assume that the 2n roots of the equation 


|Ao®?+ Bo+C|=0 
can be distributed in two sets, of n roots each, 
Sr (Cigna, wes SO ee 


such that none of the roots in S, equals any root in S,, and such that there are 
n independent right eigenvectors &,,..., E,, corresponding to the o1,...,0, and n 
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independent left eigenvectors m,...,,, corresponding to o{,...,0,. Put 
OpROyee 0 om O. £0 

(4) A,=| % % +--+ 0 aati lee Og 
\OUnOe4 6, OMLOh.8s a5 


and define two (n xn)-matrices Q, R by 


1 
(5) Ree (Eva ey ee ss 
ny 
Then 
(6) ADT BDAC = 1s (1D —A,) 1, AR IDEA) RE 


3. Proof. Multiplying (3) from the right by R, and from the left by Ly, 
we reduce the general case to the case where both R, and L, are unity matrices. 
Then each &, is the k® coordinate unity vector ex =(6,,x) Corresponding to o4, 
while 7, is also e, = (d,,,) corresponding to o;. We have then 


(Ao,* + Bo; +C)¢e,=0 (=A eye 
and, taking here the k'® row and putting 
A= (4,5),° B= (By), C= (Cys); 

(7) G0; 0, 0, Ae =O (inks =e 71) 
In the same way we have 

e,.(4o,° + Bo, +C)=0, 
or, considering the 7" column, 
(8) Gee 0.6. — 0 (ates enter aye 
From (7) and (8) we have 
(9) beg =— (C7 +0%), Cee = 44 50% 5 
this can be written as 
(10) B= As AP AA). GG = AAPA 
and these relations give (6). Our theorem is proved. 


4. Under the hypotheses of our theorem it follows from the Fredholm alter- 
native that to the roots a; in S, correspond v linearly independent right eigen- 
vectors &, ..., &, and to the a, in S, linearly independent left eigenvectors 
nw, ...,49). Further, under our hypotheses, A is not singular. 


To obtain the &”) observe that if we write 


(11) LA RPP RE AAT eT (,,), 
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to the decomposition (6) corresponds the decomposition 
(12) AE BAS CHL; (A= Ay P20 aw 
A right eigenvector &'” corresponding to a, is obtained from the equation 


13} (L 0%’ — A,) T3(I 0), — Ay) Ri" 68 =0, 
or, putting 
(13) hx (I rites Aj) Ree el) =—é£, 
(14) (Io, —A,)£=0. 
Then (14) is satisfied by £=e,=(6,,), and we obtain from (13) 
(15) ED Rela, Aja 4: 


5. From now on we use the notation A”) for the column vector formed by ° 
the »* column of A. Observe that we have then from the law of multiplication | 
of matrices 


(16) (4B) = 4 - Be), 


Further, we introduce the (7 x )-matrix 


(17) Q=( 
We have then from (15) 
Te, =T OY (a jee aor, 


‘uy (10s ae A ae tap 8) - 


Mt 7 
0, — Oy 


(18) é() — (R, Q), 


Here R,Q is not singular, and therefore (18) gives for any multiple root o, 
independent eigenvectors equal in number to the multiplicity of o,’. We see 
that the matrix 


(19) Re= BO 


consists of 7 linearly independent right eigenvectors corresponding to the eigen- 
values of S, and ordered in the same way as in S, and A,. A completely 
symmetric discussion gives us the matrix 


(20) Li =O Ls 


6. The differential equation (1) becomes, in virtue of (6) and (114), 


(21) Ls (14, — 4.) TAT — A) RYE) =C00, 
or, putting 
(22) TALS, — A) RA EQ) =a), 


(23) 13 (1 —A,)a() =£(0. 
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An integral of (23) is given, as is immediately verified by differentiation, by 


(24) a. (t) = e42* fi et Lt (tr) dz, 


and we have in the same way from (22) and (24) 
t 
E(t) = Riera“ Tals) ars 


[ee 
(25) Sa Ry eer eee err IAc (a) drat’ 


ec Cc 


Interchanging the order of integrations by Dirichlet’s rule, we have 


(26) EWS Rrew f Kye ise (sr) dx 


t 
(27) K(t) = di eh’ T ef2@ dr’. 


t 


If we consider in the matrix (27) the element in the mw‘ row and the y® 
column, we obtain for the indefinite integral, denoting by ¢,, the corresponding 
component of T, 


and this is, by (17), the corresponding element of the matrix 
e At O eet 


We have therefore from (27) and (26) 


(28) Ei) =R, f (Qe-9 — eh" 2) Ltt) dr, 
or, using (19) and (20), 
(29) & (2) =f (RRERTSULA Riest 21) (rae. 


(28) is equivalent to the formula (39) in LANcASTER’s paper [1]. 
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The Effect of a Very Strong Magnetic Cross-Field 
on Steady Motion through a Slightly Conducting Fluid: 
Three-Dimensional Case 


G'S: Se EUDFORD 


Communicated by HILDA GEIRINGER 


1. Introduction 

We continue our discussion (LUDFORD 1961) of the steady flow of a slightly 
conducting, incompressible, inviscid fluid past a fixed obstacle, when there is a ' 
very strong magnetic field H, applied perpendicular to the free stream UW). 
Previously, we considered plane flow past a cylinder of small conductivity; 
here we treat a finite insulator of general shape, whose permeability « will again 
be taken equal to that of the fluid. The present paper is, however, self-contained. 

Again the relative change in the magnetic field, due to motion-induced current, 
is small, on account of the smallness of the conductivity o, and may be neglected. — 
On the other hand, the fluid motion is no longer close to the potential flow 
which obtains for o =0, since the current, though small, experiences large pon- 
deromotive force due to the strength of MH). In the limit o+0, cH}. the 
velocity is only disturbed in the direction of the applied magnetic field and all 
quantities are constant in this direction. 

The limit equations possess a multiplicity of solutions, none of which satisfies 
the boundary conditions completely. The simplest is in general not the correct 
one. Some care is needed in deriving the correct limit equations, moreover. It 
is easily shown that if terms containing the charge density o, are neglected 
(magnetohydrodynamic approximation) the resulting limit equations lead to even 
greater indeterminateness: they ensure that the local electric field is everywhere 
zero and thereby render ineffective the boundary condition of no current flow 
across the surface of the body. Consequently, the fluid velocity across the magnetic 
lines of force, and hence the vorticity along them, is undetermined. In fact 0, 
is zero in the present approximation, as it is automatically in the plane case. 
But it is essential to know that it is accurately zero and not just negligible. 
The situation is similar to that in electrostatics, where the electric field in a 
dielectric is not determined until the volume charge is specified. In the present 
case it is sufficient to know that the fluid is uncharged upstream since this 
condition persists throughout the motion of a fluid element. 

It is at this point that our problem differs from STEWARTSON’s (1956), who 
has considered slow flow with o and H, arbitrary. Provided the ,-terms are 
neglected the problems are identical*: the same limiting values of the parameters 

* Except for the type of body considered. 
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Ry, and N [see equations (2a,b)] are involved. We then agree with his 
conclusion that the flow is not determined by the steady state equations. In his 
case it is inappropriate to include the o,-terms, since they are of higher order: 
the new parameter f? [see equation (2b)] thereby introduced is smaller still 


than R,,. In fact B?/R, = ie is the ratio of the charge relaxation time to 


. . . 0 . . . . . . 
the transit time of the fluid past the body. In our case it is large; in his it is 
small. To determine his flow he considers how the motion is set up from rest. 


This is not necessarily the reason why STEWARTSON’s steady flow is so radically 
different from ours. Another possible explanation is that the bodies are different: 
his is a perfect conductor, ours a perfect insulator. It would be instructive to 
have the solution of STEwaRTson’s problem for a dielectric sphere and of ours 
for a perfect conductor. 


At this stage we know the correct form of the limit equations, but still have to 
select the correct solution from among the multiplicity of possible ones. We are 
in fact at what was the starting point in the plane case, and can immediately 
show that certain quantities, such as the disturbance velocity perpendicular to 
the magnetic lines of force (in the y-direction), and changes along them are zero 
in the limit. To determine the solution we must, as before, compress the flow 
in the y-direction so as to magnifiy vanishingly small quantities; since all of the 
latter can now be identified, this leads to a definite transformation of variables 
and a set of “boundary layer’’ equations. 

Fortunately these equations are once more linear. They take into account 
the inertia forces in the fluid, which play a critical role at large y-distances in 
dispersing the disturbance due to the body and thereby controlling the flow 
around it. In this manner the correct solution of the previous limit equations is 
selected and the way in which the pressure becomes infinite at each finite point 
determined. 

The pressure and hence the forces on the body are of order Vo Hy (i.e. VN). 
The drag has been computed for a general ellipsoid with one axis along U, and 
another along H,. As the third axis becomes indefinitely long the drag coefficient 
tends to a value considerably less than that of the elliptic cylinder treated pre- 
viously (LUDFoRD 1961). Unfortunately, no explicit results have been obtained 
for a body which is not symmetric with respect to the x, z-plane. The underlying 
integral equation has singularities which appear to be intractable. 


2. The Equations of Steady Motion and Their Limiting Form 


If the fluid velocity is denoted by Uv, the pressure by 09 Uy p, the magnetic 
field by HH, the electric field by wU) HE, and the coordinates by (ax, ay, az) 
— where a is a length — then the flow is governed by the equations 


v-gradv =— gradp + 4, (@.E +I XH), (1a) 
J—o0,v=Ry(E+v x Hy), (1b) 
divv=0, culH=0, dvH=0, 
Bere o,=PdivE, J=curlH. 
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All quantities are dimensionless, the charge density being (H,/aU)) e,, and the : 


ters are jag Vm 
parameters 2 EAL | 0 (Alfvén number), 
is 7 
Ry, =Ujauo (magnetic Reynolds number) , (2a) 
p= U% : where c? = 1/we. 
Cc 


Fig. 1 shows the choice of axes. 
In equations (1) the terms containing g, are usually small compared to others, 
and in the so-called magnetohydrodynamic approximation they are omitted (GOLD- 
ki STEIN 1960, p. 60). However, the argu- 
| ment fails in the case of vanishingly small 
conductivity, and their omission in the 
peri ioa | present problem leads to indeterminate- 
; ness in the electric field through the 
ra charge density 0,. The latter is in fact 
zero tothe accuracy used here. In the plane 


FORD 1961) it was zero automatically, 

Take the divergence of the conduction 
equation (1b) and set R,,;=0. Then 
v- grad o,=0, so that o, is constant 


ato F-yrz on streamlines. But the incident flow is 


ietse Hota oad ace uncharged. Hence 9, =0, correct to O(1) 

in Rj,, throughout the flow. In other 

words, as far as the electromagnetic field is concerned the fluid is non-conducting, 

and for any such material the elemental charge is conserved and must be specified 

in advance. Here it is zero since each fluid element originates at infinity upstream, 
where it is uncharged. 


problem considered previously (LUD-_ 


It follows that E is a potential field which may be written in the form | 


—k-+gradqg where V29 = (3) 
then ¢ is a disturbance potential. ines H is a potential field in the limit; 
and when the permeability of the body is the same as that of the fluid this leads 
to H=j: the applied magnetic field is undisturbed. 

Now consider the momentum equation (1a) which, on substituting for J from 
equation (1b), becomes 


v- grade =— gradp+ N[x(E+v x H)+(E+vx 8H) x Hj, 
—_ Ry _ ap * Hic 
at AM Fg te U 2 I 
and x=9,/R,,. Using the results just given for E and H, we obtain for the 
component equations 


where 


ou ou Cu ou op é ie 
U ee eee eee 
Oe Oe ag ee ax +N |x Ox w) er +1), 
ov , ov , Gv ep ep 
ru ) ) = — 
ax oF aaa z Lies Taare. (4) 
0 


Fu p08 pw a By whl a) (22 wl 
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_where u,v, w are the components of the disturbance velocity v—i. When we 
| speak of a very strong applied magnetic field we mean that N is large, even though 
o is small. The equations then reduce to 


a bia seid yet os aaa 
o eas (5) 


Bee ae ee eo” 
ain dw Sea er aad 
if we allow for the possibility that # becomes large with N. To these must be 


added the continuity equation 
ou ov ow 


Ox 1 oy | Oz a (0) 
and the divergence of (1b), 
| Ox Ox Ox Ou ow 
gin pomien, 1 Oe a aie (7) 


where the right-hand side has been given its limiting value. 


3. Solution of the Limit Equations 


The following argument is similar to that used in the plane case. We assume 
that a line parallel to the y-axis cuts the surface of the body in at most two points. 


First neglect quadratic terms. Then equations (5) give 


p=Nf(x, 2), u=—fe—Qz; w=—f,+9,, 
so that (6) yields 
is v= Illex the) + 8% 2); 
here / and g are arbitrary functions of x and z alone, each of which may be dif- 
ferent in the regions vertically above and below the body. 
It is clearly not possible to ensure that conditions are undisturbed at infinity 
in the y-direction. However, it is reasonable to assume that they are only dis- 
turbed in a finite manner, 7.¢., 


fo lo =O Local Az 2" 


According to Liouville’s theorem / is constant, and if the pressure is not disturbed 
at infinity in a plane y =const. then f is identically zero. 

A similar conclusion can be drawn for gm from equation (3) once it is known 
that @ is independent of y. The latter follows from the boundary condition that 
no current flows across the surface of the body. For, according to (1b), the 
normal component of the local electric field in the fluid, 


Op op OP | ( op 
u)={0, —, 0}, 
ae Sanaa ed a ae 


must then be zero at the body; so that dq/@y is a harmonic function which van- 
ishes on the body. It is therefore identically zero. 

Finally, since « and w are now everywhere zero, equation (7) shows that x 
is independent of x and hence zero. This ensures that we have in fact the exact 
solution of equations (5), (6), and (7). 
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The complete solution is 


where for a point within the circumscribing cylinder @ the boundary condition 


on velocity at the surface y =F*(x, z) of the body (see Fig. 1) gives 


_ [Fi (%, 2), Nie (Le OMe, 


B02) =) et pE( ne. 


This violates the conditions at infinity: v does not tend to zero as y—>-+oo. 
Outside @ these conditions can be satisfied, by setting 
hy ii. 2) ihe 
However, it will turn out that this is 
not in general the correct choice (see 
end of Section 6). 

Only the verticalcomponent of ve- 


of fluid slip freely in the direction 
of the applied magnetic field on 
being displaced by the obstacle 


Fa i i, 
ED 
———_ 


(see Fig. 2). 
The electrical field inside the 
body is a potential field, and con- 


Fig. 2. Flow in plane z=const. for plausible (but incorrect) 


solution of limit equations tinuity of its tangential component 
across the surface shows that it is 
uniform. Hence E =—k throughout space and any difference in dielectric con- 


stant between the fluid and the body results in a surface charge of total amount 
zero. [This is also the solution when the fluid is replaced by a vacuum. | 

On the assumption that the pressure and electric field are continuous within 
the fluid, we have been led to the solution (8) when the body is non-conducting. 
To be sure, v is not continuous in general, but this is allowable within the frame- 
work of the theory of discontinuous solutions (voN MIsEs 1958). 

The situation is quite different when the body is, for example, a perfect 
conductor: the electric field is the same as when the fluid is replaced by a vacuum, 
so that y is no longer zero (for a sphere g =2/r3). Consequently, there is now 
motion-induced space charge x, see equation (7), arising from the component of 
vorticity in the direction of the applied magnetic field. [Of course, all this 
applies only to the linearized equations. | 


4. Compressing the y-coordinate 


As im the plane case, changes in the y-direction at any fixed point have van- 
ished in the limit N->oo and conditions at infinity are violated by v above and 


below the body. To retain these changes and relieve this violation, the y-coordi- — 


nate must be compressed as N increases and quantities which vanished previously 
must be magnified. 


@ u=w=yp=x=0, v= g(%, 2), (8) 


locity is disturbed: vertical columns ~ 


during their uniform forward motion 
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Set 
/ = U WwW D K 
—I/NY = N i —— 7) = ——_ ————— Py ipeeenterak 
y | / p V Lee u VN ) Ww VN , p VN ? % VN (9) 
in equations (4), (6), and (7); assume all new variables and derivatives are of 
order unity; and let N tend to infinity. Then 


oP od 


Wes Ox 02 U, 
oP o@ 
= (ees 
u dz | Ox W, 


and 
oK oU ow 
Ox Oz Ox (11) 


To these must be added the new form of 


equation (3), 
Dig Dea 
Axe ' az " 


The occurrence of the term ¢dv/0x in equa- 
tions (10) means that the inertia of the fluid 77 Vee) 
is now taken into account. The choice (9) Fig: 3. sonal ar et rane of equa- 
ensures that, in the new variables, the flow ne 
is the result of a balance between the inertia forces and the stresses (pressure 
and Maxwell). 

From this last equation we conclude that ®=0 in every plane Y=const. 
Then from equations (10) we find 


oP oP 
a av ___aP a EP er Geb) 
Ox Bye? Ova On? 022’ 


while equation (11) gives ¢K/0x =O so that once more the charge density K is 
zero. 


Eliminate P from equations (12b) to obtain 


C2 o2 Ov o2u 
(F: se) AS ean ame (13) 


The boundary conditions are 


Fe (x,z) for Y>+0 
—= 


on &, 14 
Pac) Stor ee eo 


where . is the area bounded by the trace of @ in the x, z-plane (Fig. 1) and 
represents the whole body in the x, Y, z-space (Fig. 3). Once v has been found 
the pressure P is uniquely determined by (12b) and then the x- and z-components 
of velocity from (12a). 

The equation (13) will be solved separately above and below the x, z-plane. 
Each partial solution must satisfy the appropriate one of the conditions (14) 
on . and lead to the same values of v and P on the remainder of the x, z-plane. 


7 
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We may restrict our attention to the half space Y>0O and suppose for the 
moment that v is prescribed everywhere on the x, z-plane. The complete solutions 
can be built up from the source solution, 7.e., the solution having 0(x) 6(z) for, 
its values in the x, z-plane. Here 6(x) is the Dirac delta function. In the plane 
case the source solution is the x-derivative of the solution having the unit step, 
function for boundary values, and the latter is a function of x?/Y? alone which 
is easily determined by Laplace transform methods. In the present case there is: 
no such solution and we must proceed directly. 


5. The Source Solution and its Properties 
Since the Fourier transform of 6(x) 6(z) with respect to x and z is 1 (LIGHTHILLI 
1958), the function 


fe.e) fo, 2) — 
vo(x,Y,2)=f fexp[nx+pz—\V—n(n?+p2)Y]dadB (n=2nia, p=2278),. 
—co —0o 
is the source solution of equation (13). For the coefficient of Y we have takent 
the root with negative real part. The root with positive real part leads to expo+ 


nentially large values as Y—-+- oo. | 
Of main interest is the corresponding pressure field | 


moh dl =f fate ae a a \— n(n? + p%) Y] da dB. 


According to ERDELYI ef al. (1954, p. 17) the B-integral can be evaluated explicitly, 
and 


B(x ¥,2)=4 f Ymer* Kol2a|al (+ YH] da. 


The remaining integral can be evaluated for Y=0, which is the case of greatest 
interest. On writing the resulting integral as the sum of Fourier sine and a: 
Fourier cosine integral, we find from ERDELYI e¢ al. (1954, pp. 49, 106) 


LAG A ive) Fw HC : f Big Dept: 
age baths 16|z|8 | 2? (2) Aueaey 7) 
5 


2 x 5 3 x2 
F( Seca 
Yale ae Cae teal te =) 


This is expressible in terms of the Legendre function of order one half (ERDELY! 
et al. 1953, p. 128); alternatively complete elliptic integrals may be used *. 


(15) 


= Aree |2E (|cos 5) — K{ cos >|)] = 7 70) (say), 


where 7, # are polar coordinates in the x, z-plane (Fig. 3). 


* This is easily verified. The only explicit reference the author has found is the 
index of tables by FLetcHEr, MittER, & ROSENHEAD (1946, p. 238). 
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Fig. 4 shows a graph of A(#) for 0<%<_q; the function is even in #. For 
# =7/2, t.e. on the z-axis, 


P= = 0.0761. 
| ery = 97 
As 3-0 (positive x-axis) 0.16 — Fr 
Ole 
Pr ae log sin#?, 0.08 
OA fad He 
| , 
while for =a (negative py % 
X-axis) Bear 
, —008 
Y Sega sects A 
! = 016 
The logarithmic infin- — 020 
ey*in 7 forf=0is curi- ~7# t 
ar Bo : D0 20 30 0 50 60 1 80 3 100 0 120 130 140 150 160 110 180 
ous, since in the two-dimen- 0/7 degrees 


sional case of a line source 
along x=O there is no 
downstream influence, 7.¢., the pressure is zero for —2/2<0<a/2 (LUDFoRD 
1961). However, the present solution checks with this previous result: when 
integrated over a line of sources, JJ gives the correct two-dimensional pressure 
(Appendix 1). 


Fig. 4. The pressure function A(#) associated with a unit source [equation (16) ] 


6. Determination of the Complete Solution 
The solution of (13) taking on the values V(x, z) as Y>+0 is 


Dita, 2) Se (frogt (x—&, Y,z—n) V(é,n) d& dn, 


where the integration is to be taken over all values of € and 7. There is little to 
be gained by writing down formulas for the corresponding P, U, and W. All we 
need is the value of P as Y->+0, and this follows immediately from (16). To 
satisfy the first of the conditions (14) we set 


Vilas FS (xe2)erry ine? 


there remains the determination of V on the remainder of the x, z-plane. 
As Y->-+0 the pressure is given by 


P(x, +0, 2) ee 2—n) Fy (8, ) d& ant JJ Mae, 2—n) V(E,n) dE dn, (17) 


where C(./) is the remainder of the x, z-plane. For the corresponding solution 
in the half space Y<0, the pressure as Y->—0 is given by the same formula 
with /7 replaced by —/I and F; by Fy. Use of the same function V(x, z) in C(.7) 
for the two cases ensures that v takes on the same values as Y-++-0 there. To 
ensure that the pressure [and hence the remaining velocity components] are 
continuous across Y —0 outside ., we must take 

Sf I(x — & z—m) V(E,n) 4& dn 


C(f) (18) 
ee ee — §,2—n) [Fi (6,0) + Fy (8,9) do dy 


* This arises from K(k), which is logarithmically infinite for k = 1. 


250 G. S. S. LUDFORD: 


for all values of (x, z) in C(/). Having determined V(&, 7), we may calculate the 
pressure on the upper and lower surfaces of the body from (17) and its analog ' 


For a body symmetric about the x, z-plane f° = — F; and the solution of (18) is; 
V4, zb= 0) anCiZ); 


which is otherwise obvious. The integral equation has so far not proved tractable 
for any region ..7 when the body is not symmetric. In any event it is clear that V 
will be non-zero. 


7. The Drag on an Ellipsoid 
The drag experienced by the body is 


D=0,U2%a2\N Jf [Fi (x, 2) P(x, + 0,2) +E (x, 2) P(x, — 0, 2)]dxdz, 
a 
and for a body which is symmetric about the x, z-plane this becomes 
D =29,U2a?|/N ff dxdz ff d&dnF; (x, 2) I(x —&,2z—n) F7 (En). 


For an ellipsoid F*(x, z) =k\/1— x?— «222, where a, ka, a/x are the semi-axes,, 
and D =20,U?a?\/N k?I/2a, where 


Tg yee i dxdz [| dé dy aL Col I te 
Tt A V(4— #2 — 2 22) (1 —E2— 0? 7?) 
x? to2 22S E24 G2 y2<1 


NT is the drag coefficient of the spheroid k =1 based on the area za. (a/x) of . 


—— lylinder The integral has been evaluated 
numerically (using an eight digit, — 


ce floating decimal scheme) for «% = 
32 0(0.2) 1.0(0.5) 2.5 and a table is 
80 given below. The convergence with 
| decreasing mesh size was slow for 
oe small values of « and the results 
| | had to be extrapolated*. A graph 
as of J versus « is shown in Fig. 5, 
22 The value J(0) =3.36, correspond- 


20h > ing to an infinitely long transverse 


20 ens : 
Shiri oie. Alby % axis, 1s considerably less than that 
ig. 5. Values of the i al (4 termining dr ipsoid ‘ ; 
ntegral (19) determining drag on an ellipsoid of a circular cylinder, Bah, 4)/ 
182 = 3.66 (LupFoRD 1961). These values compare with I(1) = 2.76 for a 
sphere: 


a=0.0 0.2 0.4 0.6 0.8 1.0 1.5 26 25 
Ie) =3.37 3.30 3.46 4.02 2.89 276 2.48 DOS 2.10 


| 
* See Appendix 2 for the reduction of I(x) to a form suitable for numerical | 
computation. | 
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Appendix 1 


| The pressure exerted at the origin by a line of sources of unit density along 
X% =X) >0 is, according to equation (16), given by 


$d, 


P(X) rhe Xp, 2) dz ~ aa] 70) |sec 0 
where z =— x, tan #. 


Now A(#) may be written as the sum of even and odd functions of cos #, 


Pp dae 3 ge a a 2 _ 2cos$ Cy ce 26)].* 
@) 16 rp Fla 7?) 008°9) I) Fp, gg 0088). 
Also, with cos ? =— «?, 


fF, eer = >i cos?) |sec d|! dd 


1 
ae 3: ac et 6 ee fees BP 
Safe ea Veeet Bs 4, Te 
0 


See 18 3% 4 
fF. oy 808°) [cos dad 


i 
1 LP CES aay ee er Te eRe IED 
=5 fF. 4.2 :0)« (14 —a)~*da Tales) 
0 


These last results may be obtained by use of the integral representation (ERDELYI 
et al., 1953, p. 59) 


1 


EGROSC: h) ORM P24 — 2-1 (1 — ta) "dE, 


0 
Re (c) > Re(b)>0, and interchange of orders of integration. In the former case 
the hypergeometric function must first be written as 


aL 
1— af F(Z, 8, 3:05) ds, 
0 


1.é., as the integral of its derivative. 
It follows that 


eo alt 1 1 
P(X) Vx aya | Tal ines for x) > 0. (21) 
Similarly n/2 
P(x») = ee (sec#)#d9=0 for x%<0, 
i aA) 


since the substitution cos? =«* leads to the same «-integrals (20) and hence to 
(21) with + changed to — in the bracket. 


* This is the result of applying the transformation 
TAG, (oy, (03) = (Gl) 7 JO(7, Od, (2 41) 


(ERDELYI et al. 1953, p. 64) to the hypergeometric functions in (15). 
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Appendix 2 
The integrand in (19) is singular not only on the circumferences of the ellipses 
but also for z= when x=é (and especially when x =). In order to transpose , 
the latter to the boundary of the region of integration, set (Fig. 6a) 
x=A+rcos8#, &=A—rcos#, 
z=p+rsind, y=p~—rsind. 


(a) (b) 
Fig. 6. (a) Change of coordinates x, z, §, 7 to A, u, 7, 9. (b) The integration region A(r, #) for A, u 


Then (Fig. 6b) for fixed 7, # the point (A, w) ranges over the intersection, A(r, #), 
of the two ellipses obtained by displacing the original through (+rcos#?, +,rsin#). 
Thus 


I= 1s [20 a0 all oe didu 
R(A, u37, B) : 


where 
R(A, “3 7, 9) 
=[1—(A+rcos#)?— «(w+ rsin 9)?] [1 — (A —rcos 8)? — a2 (u — rsin d)?] 


and /(#) is defined below. By using new coordinates X=A cos 0 +02u sin 0, 


Y=—aAsin#-+-«u cos # and integrating over Y, the integral over A may, be 
reduced to 


= [ (BX, 8) K(c) — C(X;7, 9) E(c)] dX 


where 


f=co?d+oa2sin?d, g—o2sin2d, o = | eri oar 


_ [f—2g) X*+2rfgX+f(e—1?7?)] 
PYi—(X—r fy? a 


Cond pace eerie: 


and K and E are complete elliptic integrals of the first and second kinds respec- | | 
tively. 
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The limits of integration are fixed for the variables #, s, ¢ where ry =/~s and 
_X=(1—s) ft. The logarithmic and/or algebraic singularities in the resulting 
integrand of J at #=0, s=0, #=0 and t=1 can be removed by setting 9 =¢?, 
f—o', ¢=sin2y. 

For «=0 there is a further singularity (arising from f=0 at ? =2/2). How- 
ever, the #-integration can then be carried out analytically using the results 
given in Appendix 1, and the remaining o- and y-integrations are regular. 


| The numerical work was carried out on the IBM 7070 at the Brown University 

Computing Laboratory. Thanks are due to Mr. Gorpon T. Witmot who carried out 
the programming. 
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A Uniqueness Theorem in Magnetohydrodynamics 
R. H. DyER & D. E. EDMUNDS 


Communicated by D. GRAFFI 


Introduction 

Consider the flow of an incompressible, viscous, electrically conducting fluid: 
around a finite solid body B, which is assumed to have the same magnetic: 
permeability as the fluid. It is supposed that there are no other boundariess 
in the fluid, which extends to infinity in three-space. We show that subjectt 
to certain conditions on the velocity and magnetic fields v and H and on theirt 
first spatial derivatives, the fluid motion and the electromagnetic quantities are: 
uniquely determined by the motion of B and the distributions of v and H att 
some initial instant, together with the behaviour of the pressure p at large: 
distances from B. 

The method used is an extension of that used recently by GRAFFI in a series = 
of papers [/, 2, 3] dealing with the corresponding uniqueness theorem for the : 
non-magnetic case, both incompressible and compressible fluids being considered. . 
An important feature of these papers is that uniqueness is proved without : 
assumptions about the convergence at infinity of the velocity, such assumptions + 
being common in earlier papers on this subject (e.g. [4]). Uniqueness theorems ; 
for magnetohydrodynamic flows have been proved by NARDINI [5] and KANWAL . 
[6] under various hypotheses; FERRARI [7], following the method of GRAFFI, 
has established uniqueness in the magnetohydrodynamic case under conditions : 
of convergence at infinity weaker than those required by NARDINI. 

The present result is similar to that of [7] as far as conditions at infinity’ 
are concerned, but differs in the magnetic boundary conditions used on the 
surface of the solid 6. FERRARI assumes that the tangential component of the 
magnetic field is specified on the body surface at all times; here all we require 
is that the magnetic field and the tangential component of the electric field be 
continuous across the body surface. 


Equations of motion 


The equations governing the motion of the fluid are, under the usual as- 
sumptions of magnetohydrodynamics [8], 


Ow Mu £: 
Om rael o(v- grad) » — _ (H- grad) H = — grad @ + ovV2v+ oX, 1 


47 ( 

OH ; ad) H 2 
3p + (v- grad) H — (H- grad) v = V? H, (2 
curlH = 4z0(E+yuv~x HH), (3 
div u=0, (4 
( 


) 
) 
) 
) 
divi = 6, 5) 
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where v is the fluid velocity, # the pressure, @ the density, » the kinematical 
viscosity, H the magnetic field, E the electric field, ~ the permeability, o the 
fluid conductivity, 7=(42 0) the magnetic viscosity of the fluid, X the 
externally applied force, and ©=f+H?/82. 

Let o be the conductivity of the body B, and 7=(426)7 its magnetic 
viscosity. The equations satisfied by the magnetic field in B will be (2), (3) 
and (5), with o replaced by 6 and 7 by 7%; in these equations v will denote the 
velocity of a typical point of B. 

B is assumed to be a bounded closed region of three-space, having a piece- 
wise smooth surface S. The boundary conditions which must be satisfied on S are: 


v, H, and the tangential component of E are continuous across S. (6) 


Let D be the set of points of three-space exterior to B, and let T be the 
closed time interval T= (0, ¢)), where ¢t) (>0) is arbitrary but fixed. To make 
the initial value problem precise we shall impose the following conditions on 
the various fields defined interior and exterior to B; in these conditions x; 
(¢=1, 2,3) will denote rectangular Cartesian coordinates measured from some 
origin 0 fixed in space, and 7?=x;x,; (using the tensor summation convention). 

(1) The velocity components v; and their first derivatives with respect to 
the x; and to the time ¢ are continuous bounded functions of x; and ¢ for all 
points in the product set DT. The second order space derivatives are con- 
tinuous for all points in Dx T. 

(u) The components H; of the magnetic field and their first derivatives with 
respect to the space coordinates and the time are continuous bounded functions 
of x, and ¢ for all points in D xT and for all points in Bx T. (H; is continuous 
across S by (6), but there may be discontinuities across S in 6H,/éx,;.) The 
second order space derivatives are continuous for all points in D xT and for 
all points in BxT. 

(iii) The pressure # is continuous and has continuous first derivatives with 
respect to the x,;in DT. At infinity p converges to a given value fy such that, 
uniformly in ¢, for all f€T, p= Pp, er) 


for large 1. 

(iv) The v, are specified throughout D for t=0, and the H; are specified 
throughout DUB for t=0. 

(v) The v, are specified on S (and hence, of course, at every point of B) 
for every ¢20. 

(vi) The external force components X; are specified as functions of x; and ¢. 

Put simply, we envisage a situation in which the body B moves with pre- 
scribed and finite velocity in an unbounded fluid in which the pressure approaches 
a specified value fy at large distances from B; in addition the values of v and 
H at some initial instant are prescribed throughout all space. The content of 
this paper is summarised in the following 

Theorem. There can be at most one solution of equations (1)—(5) (and of the 
corresponding equations which hold inside B) satisfying the conditions (6) and (1) — (vi). 

Proof. Suppose there are two solutions v;, H;, pf and v,+u;, H; +h, p+f' 
which satisfy equations (1)—(5) and conditions (6), (i)—(vi). Then at points 
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in the fluid, 


3 e F) oh; 
Ou; ta) ’ oy, b :: a 
at ie, OX, (vj + j) + OX, 4n| * Ox, bit y Py lioreda OX, 
00)’ OU; 
gg Ear ton, 
, Gu; ah 
oh 0 iy Vappsdec sees i 
a r Uy ax (A; + hj) + % é hy é (0; + %) R Ox OX, OXp (8) 
oui 9, (9) 
OX; 
eh 9, (10) 
Ox; 
mancre @! =p’ + (u/82) (2H; hj +h; hj). 


For points lying inside B the corresponding equations are 


a “hj 8! 
at a er My OX, Ox, OX, (8) 
Ue eG (10’) 

Ox; 


V=(V,) being the (prescribed) velocity of a typical point of B. 
Let wp be the closed ball of centre 0, radius R, and of surface og. We choose 


R so large that for all ¢€T, BCwpr. Let Wpr=WpoCB, where cB denotes the » 


exterior of B. Then on multiplying (7) by u; and summing, integration of the 
resulting expression throughout wp gives, after using Green’s theorem and the 
boundary condition u;=0 on S, 


Ou; Ou; 


1 
il le ay (jm) Bee Ox, OXp 


7) eas ) oh; 
—— of jum om (v; + u,) ie uy {hy Bx (H; + h,) Eager} 


AW p 


dwp (11) 


Ou; ast q 1 : g 7 
QV Uj a, — @' Uj ;| donp— > O | Uj Uj Uy My AOR. 
OR OR 
Here n=(n,) is the normal drawn outward from wp. Similarly, if (8) and (8’) 


are multiplied by 4; and summed, then integrated throughout wp and B respec- 
tively, addition of the resulting expression gives 


1 2 hy eh; . oh; ah; 
[la ae umd tat olde t fF aed hag 
B 


dB 


OX, Ox, 
maf he moligers | (1-4) mas — 


Fs) Ou; 
site fae (H; + h,) +54 cae (v; + 4) — Hy "| dWp— 


fa (4st nay ot) aB. 
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, In obtaining this the boundary condition (6) on S has been used; A(f) is the 
_ discontinuity in f on crossing S. Combining (11) and (12), 


) lle (5 Qu, u; + i hi) +e ue 4 nM oh; oh; ioe : 


ot OX, OX, 4m OX, OX, 
OR - 
a ~ dh; ah; 
| + fae lara) tae, | 2 
/ 
foal 1 
| ect te acres 
Wp OR 
mh Ou; 
— [|e G aaa Oe pepe" dop+ 
oR 
‘ (13) 
lu B) oh; 
tHe fifties +h) + Ha AN 
OR 
Ou; F) oh; 
+h, {he 5 (v; + 4,) + A, a. ox, (H; + h,) — v, ve \| don — 
lu = Ny) 
= filte - hy ‘\aB 4 
st fag mdort ae fWAlrg FA) ma S. 
Now 


oh; ol 
[An gi)mas = fh Alqeuthxnjds + fh; 4 (n aa) mS. 
S eS) S 


Denoting the electric fields in the two postulated possible flows by E and E-+ E’, 
we have, in the fluid, 


E=—wvxH-+ (420) 1 curl 
E+F=—wu(v+u)x(H+h) + (400)7 curl (H+ h) 
so that E' = (420) curlh — uw{ux(H+h) 4+ (vxh)}. 


For points inside B a similar expression for E’ exists, with o replaced by o and 
wv by V. At a solid-fluid interface the tangential components of the electric field 
must be continuous, and this provides the condition 


O=nXxA(E’) =unxAlycurlh —ux(H-+h) — (vxh)]. 
Since » and A are continuous across S, and w=0 on S, this reduces to 


Afncurlh] xn=0. 


fr Al oi) maS = fh, A (n*) mas 
Ss S 


== fh Fem do nf [e oe (curl h)?) dwp+ 
7 


Hence 


é Ox, Ox, 
OR Op 
_ Pf ah; ah; : 
Li | ae dae ~ (curl h) |aB 
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using Green’s theorem. Equation (13) thus becomes 


OU; OU; 


Mla a ll zB 
| E Hei t+ aa hy hy) + ov Oxp 


OX; 


Se 1B 
dong +o [53 a © (hj h,) d 
wR B 

i 0 1 
=— 0 | uj Up (v; lL U;) dW p — a 0 f U; U; Up Ny dp — 


or oR 
= | (7 U;N; — OV U; a ny) dop+ 
+t i (04; ie — 04) so (Hy + hy) + ayy a 4 «)|dox+ (14) 
Oe De hj hj ,| do + 4 f hie sg aB + 
° B 
2 4( Pm don 
ae (curl h)?dwp— us f (cul heaB. 
we ie 


Integrating (14) with respect to ¢ from 0 to ¢,, and then again with respect 


to t, from 0 to b (St,), where 6 is a strictly positive constant to be specified 
later, there follows 


b 
seach pietlt oa atalaaes 
b 
au; 8 
serfauferf a ae den wr fan f & 
ga 
b ty 
nu 7] 
ue fas TE NERCR ES ii ea a h)2dB 
~ma ft fat fy ya (Y + U;) dp — sofas [4 [meinen 


j4;d B+ 


(15) 
fan faef U;N; — OVU; i my) dog + 
0 OR 
a a 
iG O 0 
| at, [ at [ \(u, fy — ij ty) oo (Ht Wy) + i Py ge (j+ 4) |dog+ 
0 0 oR 


b t, 
aaa Jl u; h; — 5 Hi hj) don + 


+f feferl fi hg mae Ue fan fa f( oe) tude 
k 


use being eas of ie initial condition tiv): 
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Let 
é 
(a) , 
M =sup| fag Hi +h), -M’ =sup|H,| 


the suprema being taken over all points in D x T and over all 7, k; let 


the supremum being taken over all points in B x T and over all7, k. Then using 
Cauchy’s inequality repeatedly, we find that for points in D x T, 


7) 
U; OF rm (v; + 4,;)| S3N 4,4, 


|u; u; 0, 2, | S3.N' uu, 


Ou; 3 1 Ou; Ou; 
ey pS ae 
1 Ox, UN e2 eR D (iG, (O73 
hg AH eh,\| SMa tu, hap h,) 
Ce ie |e TN Ma 


S3Nh fh; 


|H, u; h; n,| SS M'(u; u; +h; hj) 
|; hj 0, 2,| S3N’ h; h; 


(a) h; ra) hy, 
hs ae oe) 


< hy h; + (curl h)? 


| 4, Ay, Uj 2;| S3 ON; h; 
while for points in Bx7, 


OV; A 
h; hy, ae =< 3N h; h;. 


Finally we need an estimate for f p’u,n;dog. Using Schwarz’s inequality, 
oR 


Sp undo <| Spr dog: fujujdop 
OR OR OR 


ind since ~’=O(r) by condition (ii), it follows that there exists a strictly 
ositive constant c such that 


1 
2 


Se 


Sp’ un; dor Ju; u; dep 
OR oR 


ovided R is large enough. 
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Using all these inequalities in (15), we find 


b 
M 
fe) fan famton’ 


0<G(R) <34(Net 


0 Op 
b 
3 ‘ M’ ; oe ou 
+30(N etort VE [dt [ uudort Fe fas saith 
0 OR 
+34 b(M+N) fap hy ding + 34 b(2M'+N'+n+0) [esfandoe (16) 
OR 


ee pee fas fn 4 dB + My fanfaft (curl h)?dop+ 
0 
eff [omen 
0 0 OR 


where 


b 
n= fap al (omy + Ho hyh,) doet f dt, [ fh hyaB + 
0 B 


cer feafe Be 9 (17¥ 
OX, 


OX, 


b ty 
+28 f a anaes dog+ i [ dt, f at [ (curlh)? dB. 
0 Gy aie? 


0 0 Op 


By Schwarz’s inequality, 
3 b ty ty 4 
= fas, {Sach { fae f u,u; doe} 
0 0 0 OR 
b b 4 
so|{ fai} { fat fujujdox} 
0 0 OR 


b ty : 
fat, fat| Su; U; dop 
0 0 OR 


r 3 

ss al f dt Susu; dog e: 
0 6 
Putting F 


mm, = max [oN + SME 9G( SAN), oN”| 
27 0 


and choosing b=max {to/n : to/n<(2m,)+, n integral}, (17) becomes 


b 
1 3 t M’ 
G(R) S>G(R) +3 (N otvet SA) f dt, [ujujdor+ 


0 o 
b t = 


1 b 
1 Sy Cu; Ou; 3 ; : 
+ yen fade, fae fo oq, don + 34 (2M +N ++) [ dt, [hjh;don+ 


0 OR 
b 1 b L 
+ 4h fat, fat | cur)don + 28 | ttf mm dom 
0 0 oR 0 oR 
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On setting 


m= 2max[35(N’-+y+2°4), 5, 350M’ +N’ +y+Q) 
we find 


G(R) SmG'(R) +1[G'(R)}* (18) 
Where / = 2c b3(2/0)3, G’(R) =dG/dR. 


We propose next to show that G(R) =0; the argument is exactly the same 
as that given by GrarrFi [2], but we reproduce it here for the sake of com- 
pleteness. 


In view of conditions (i) and (ii), 
G(R) 0 (R®) (19) 


for large R. Suppose there exists an Ry such that G(R,)-0. Since G(R) is 
always non-negative, this implies that G(R,)>0. Because G’(R) is also non- 
negative, G(R) is monotonic increasing, and hence 


0<G(Ro) S$ G(R) SmG'(R) +1[G'(R)}* (20) 


for all R>R,. It follows that there is a constant «> 0 such that [G’(R)]}*=« 
for all R>R,. Then (20) may be written 


G(R) <| Te G(R) S(m+ 1)G"R) 
that is, 
aa alnt =) =5>0 
and so 


for all R>R,. This contradicts (19), and hence G(R) =0. We deduce that 
u;—0 throughout D x (0, 6), and h;=0 throughout (DUB) x (0, 8). 


Now integrate (14) with respect to ¢ from 5 to b+4,, and then with respect 
to ¢, from 0 to b. The previous arguments may then be repeated to show that 
u,;—O0 throughout D x (6, 2b), and that 4;=0 throughout (DUB) x(b, 20). Pro- 
ceeding in this way we eventually cover the whole of the arbitrary interval 
T = (0, ¢)) in steps of length 0, so that we have finally 


“u=0 m DxXT, h=0 m (DvB)xT. 


It follows easily from (7) and (iii) that #’=0 in DXT, and hence the solution 
is unique for all ¢€ 7. But T is arbitrary, and hence the theorem is proved. 


We are most grateful to Professor GRAFFI for a suggestion which led to a shortening 
of the proof, and for drawing our attention to references [2, 3, 7]. 

One of us (R.H.D.) is indebted to the Department of Scientific and Industrial 
Research for a maintenance grant. 

Arch, Rational Mech. Anal., Vol. 8 18 


262 R.H.DyvEerR& D.E.Epmunps: Uniqueness of magnetohydrodynamic flows 


References 


(Z| Grarri, D.: Sur un théoréme d’unicité pour le mouvement d’un fluide visqueux : 

; dans un domaine illimité. Comptes Rendus 249 (18), 1741—1743 (1959). 

[2] Grarrt, D.: Sul teorema di unicita nella dinamica dei fluidi. Annali di Matematica | 

(LV) 50, 379—388 (1960). 

[3] GraFFt, D.: Sul teorema di unicita per le equazioni del moto dei fluidi compressi- 

bili in un dominio illimitato. Atti della Accademia delle Scienze dell’Istituto | 

di Bologna (XI), 7, 1—8 (1960). 

[4] FoA, E.: Sull’ impiego dell’ analisi dimensionale nello studio del moto turbolento. 

L’Industria 43, 426—429 (1929). 

[5] Narpint, R.: Due teoremi di unicita nella teoria delle onde magneto-idrodina- 

miche. Rend. Sem. Mat. Univ. Padova 21, 303—315 (1952). 

[6] Kanwar, R. P.: Uniqueness of magnetohydrodynamic flows. Arch. Rational 
Mech. Anal. 4, 335—340 (1960). 

[7] Ferrart, I.: Su un teorema di unicita per le equazioni dell’ idromagnetismo. 
Atti del Seminario Matematico e Fisico dell’ Universita di Modena 9, 205 —217 
(1960). 

[8| Cowxtne, T. G.: Magnetohydrodynamics. New York: Interscience Publishers 

1957. 


Mathematics Department 
University College 
Cardiff, Wales 


(Received May 27, 1961) 


EDITORIAL BOARD 


H. BECKERT 
Mathematisches Institut 
der Universitat Leipzig 


R. BERKER 
Mihiirdar Caddesi No. 104 
Kadikoy, Istanbul 


L. CESARI 
University of Michigan 
Ann Arbor, Michigan 


EACOLUATZ 
Institut fiir Angewandte Mathematik 
Universitat Hamburg 


A. ERDELYI 
- California Institute of Technology 
Pasadena, California 


JE LE RICKSEN 
The Johns Hopkins University 
Baltimore, Maryland 


G. FICHERA 
Istituto Matematico 
Universita di Roma 


R. FINN 
Stanford University 
California 


HILDA GEIRINGER 
Harvard University 
Cambridge, Massachusetts 


A H. GORTLER 
Institut fiir Angewandte Mathematik 
Universitat Freiburg i. Br. 


D. GRAFFI 


= Istituto Matematico ,,Salvatore Pincherle“‘ 


Universita di Bologna 


A. E. GREEN 
King’s College 
Newcastle-upon-Tyne 


J. HADAMARD 
Institut de France 
Paris 


L. HORMANDER 
Department of Mathematics 
_ University of Stockholm 


C. TACOB 
Institut de Mathématique 
Académie de la République Populaire 
Roumaine, Bucarest 


M. KAC 
Rockefeller Institute 
New York 21, New York 


tT KATO 
University of Tokyo 
Bunkyo-ku, Tokyo 


E. LEIMANIS 
University of British Columbia 
Vancouver 


@ CzEIN 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 


W. MAGNUS 
Institute of Mathematical Sciences 
New York University 
New York City 


G. CeMcVirE te 
University of Illinois Observatory 
Urbana, Illinois 


J. MEIXNER 
Institut fiir Theoretische Physik 
Technische Hochschule Aachen 


C. MULLER 
Institut fiir Reine u. Angewandte 
Mathematik 
Technische Hochschule Aachen 


W. NOLL 
Carnegie Institute of Technology 
Pittsburgh, Pennsylvania 


A. OSTROWSKI 
Certenago-Montagnola 
Ticino 


RS. RIVELN 
Division of Applied Mathematics 
Brown University 
Providence, Rhode Island 


M. M. SCHIFFER 
Stanford University 
California 


J. SERRIN 
. Institute of Technology 
University of Minnesota 
Minneapolis, Minnesota 


E. STERNBERG 
Division of Applied Mathematics 
Brown University 
Providence, Rhode Island 


R. A. TOUPIN 
Naval Research Laboratory 
Washington 25, D.C. 


CURUE SDE 
The Johns Hopkins University 
Baltimore, Maryland 


H. VILLAT 
47, bd. A. Blanqui 
Paris XIII 


CONTENTS 


Jenkins, H. B., On Two-Dimensional Variational Problems in Para- 
metric Form . 


BERGMAN, St., Some Properties of a Harmonic Function of Three 
Variables given by its Series Development . 


GILBERT, R. P., A Note on Harmonic Functions in (P+2) Variables . 


SEGEL, L. A., Application of Conformal Mapping to Boundary Pertur- 
bation Problems for the Membrane Equation 


OstrowskI, A. M., On Lancaster’s Decomposition of a Matrix Differen- 
tial Operator 


LupForD, G. S. S., The Effect of a Very Strong Magnetic Cross-Field 
on Steady Motion through a Slightly Conducting Fluid: Three- 
Dimensional Case Bie Rao 


Dyer, R. H., & D. E. EpMunps, A Uniqueness Theorem in Magneto- 
hydrodynamics at at a. See 


Druck der Universititsdruckerei H. Stiirtz AG., Wiirzburg 
Printed in Germany 


» £Ot 


- 207 


223 


- 228 


. 238 


> 24 


. 254 


~ 
eee Ar 


ee ae eee fee ee 


ee eS ee 


po er ee 


